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CHAPTER I 


Introduction 

1.1 Objective and Motivation 

In this thesis, the electromagnetic scattering from a cone frustum is investi- 
gated using the high-frequency techniques of the geometrical theory of diffraction. 
The emphasis in this work is on the development of a solution which, besides being 
accurate, is easy to apply and fast to calculate. In other words, we want a solu- 
tion suitable for engineering applications. Even though the scattering from basic 
shapes, like the cone frustum, is fundamental, there is much to be gained from its 
study. 

The solution to such a problem can give insight into the development of solu- 
tions to other scattering problems. There are many subtleties that need to be con- 
sidered when the solution to the scattered fields of such a geometry is attempted; 
the experience gained may be applied to similar situations in other geometries. 
In fact, the topic of this thesis was inspired by the recent work of Kuei-Chien 
Chiang (l) who investigated the scattering from a cylinder. 

From a more practical standpoint, the study of high-frequency scattering from 
basic geometrical shapes has applications in the modeling of complex scatterers. 
Due to the highly local nature of high-frequency scattering, the solution to complex 
scatterers may be modeled approximately as a superposition of the scattered fields 
from simpler shapes — such as cone frustums, cylinders, ellipsoids, plates, and oth- 
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ers. A fast and accurate solution for these basic geometries may lead to analytical 
computer-generated scattering patterns for aircraft, ships, reflector antennas, or 
even a space station. This can result in the computer design of scatterers to meet 
various engineering specifications. 

1.2 The Problem 

This is an electromagnetic scattering problem. The objective is to find the 
scattered electric and magnetic fields for a known incident illumination of the scat- 
terer. The scattered fields consist of those fields that are reflected and diffracted 
from the scatterer. For this problem let us specify the following conditions: 

• Shape — elliptical right cone frustum (see Figure 1) 

• Material — perfect conductor in free space 

• Size — greater than a wavelength (high frequencies) 

• Illumination — a plane wave polarized in either of the orthogonal polarizations 
(vertical or horizontal) 

• Observation— bistatic, in the far field with matched polarization 

The variables in this problem are 

• Frequency 

• Polarization 

• Positions of source and receiver 

• Cone frustum dimensions 
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Figure 1: An elliptical right cone frustum 
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1.3 Solution Strategy 


To solve our problem, we depend a great deal on the work which has been 
done in the past. We mentioned previously that a techinique called the geomet- 
rical theory of diffraction (GTD) would be used. It was originally developed by 
Keller [2], and it has been improved by many others which include the work of 
Kouyoumjian and Pathak [3]. The improved version of GTD, which corrects the 
areas where Keller’s original theory had failed, is usually referred to as the uniform 
geometrical theory of diffraction (UTD). A brief background on UTD is provided 
in Chapter 2. 

In the past, GTD has been applied to the scattering from a circular cone frus- 
tum by Bectel and Ross [4,5,6] and Ryan and Peters [7|. The GTD solution consists 
of a superposition of high-frequency semi-infinite wedge solutions associated with 
each local scattering center at the edges. The total scattered field is obtained by 
summing the diffracted fields from each scattering center. (Chapter 3 provides a 
more detailed explaination.) By considering only first-order edge diffraction, their 
results show good agreement with experimental measurements (in the region where 
GTD is valid). 

The only drawback to the GTD solution for a cone frustum is that the edge 
diffraction coefficients lack important information about the surface curvature and 
about the presence of “rim caustics”. The information dealing with the surface 
is necessary to obtain the proper reflected field. The “rim caustics” are places 
where the diffracted field contributions come from the entire rim rather than local 
scattering centers. And to get the correct edge diffracted field, the behavior near 
rim caustics must be known. 

Due to this missing information, the original GTD solution fails in the rim- 
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caustic regions and in the region associated with specular reflection (a reflection 
caustic region) from the curved surface. Bectel and Ross had avoided these regions 
by switching to the physical optics solution for both the specular region of the 
curved surface and in the rim caustic regions associated with specular scattering 
off the endcaps. Ryan and Peters had also used physical optics for the specular 
region of the curved surface, but for the rim-caustic regions the scattered fields 
were evaluated by a technique called the “equivalent current concept” . This worked 
fairly well, however, it would be more convienient if we could somehow preserve 
the GTD edge-diffraction form of the solution in all regions — without having to 
switch from one form to another. More recently, Chu [8], based on similar work by 
Ryan [9], has modified the “caustic distance” in the GTD solution to provide the 
physical optics result in the specular direction for the curved side of the frustum. 
While sufficient for cone frustums with small cone angles, it does not provide 
a uniform solution. A uniform solution is one that provides the correct result 
uniformly (from the regions where the GTD edge diffraction solution is valid) 
through the regions of failure. 

In this work, we will try to retain the GTD form of the solution, and yet have 
a uniform solution. In our approach, we will develop “transition functions that 
multiply terms in the GTD diffraction coefficients (thereby, generating a modified 
UTD diffraction coefficient). These transition functions will come into effect at the 
caustic regions to correct the edge-diffraction solution; outside the caustic regions 
they will have no effect. This approach has certain advantages in a computer 
code. The transition functions will eliminate the need to decide when to switch 
from one solution to another. The simplicity of the GTD edge diffraction solution 
is preserved with only a slight modification to the diffraction coefficients. 

In the following chapters we will show how the transition functions can be 
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obtained. We will show that by using “equivalent line currents” on the curved 
surface of the cone frustum, a transition function to correct the failure in the 
specular region (of the curved surface) can be developed. In a similar way, the use 
of equivalent edge currents on the rim will lead to transition functions to correct 
the failure in the rim caustic regions. 

For the cylinder, Chiang [l] had obtained a uniform solution through the rim 
caustic regions. Even though there are some differences for the cone frustum, her 
work will be applied (in addition to the rim caustic correction of Ryan and Peters). 

Bcisically, the approach used in obtaining our solution is as follows: 

1. The theory is developed. 

2. Using the theory, calculated results are generated for sample cases. 

3. The calculated results are compared with measurements or other known so- 
lutions. 

4. If the match is good enough then this will be considered the approximate 
solution. Otherwise, the theory will be refined, and the procedure repeated. 

This approach is illustrated as a block diagram in Figure 2. It outlines the contents 
of this report. Information about each block in the figure will be provided in the 
following chapters. Chapter 3 discusses the diffraction-point method solution, 
which provides the basic form of the solution. Chapters 4 and 5 deal with the 
corrections to the diffraction-point method solution and the development of the 
transition functions. A comparison with experimental measurements is made in 
Chapter 6, and with the moment method solution in Chapter 7. 
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Problem Statement 
and Constraints 



Figure 2: A solution strategy for the cone frustum scattering problem using UTD. 
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1.4 Notation, Abbreviations, and Symbols 


This section contains some information on the notation, abbreviations, and 
symbols which may be encountered in this report. It is provided as a help to the 
reader. 

1.4.1 Notation 

In this work, the time variation will be assumed, and it will be suppressed. 
Numbers in parentheses ( ) refer to equations. Numbers in brackets [ ] refer to 
references. 

For quantities superscripted with a ± or the top (bottom) sign corresponds 
to the top (bottom) sign of similarly subscripted quantities in the expression. For 
e.xpressions containing quantities with subscripts like, the first (second) sub- 
script corresponds to the first (second) subscripts of other quantities in the expres- 
sion. Also for expressions like: Rg ^ i the first (second) subscript corresponds 

to the top (bottom) sign. A quantity with a named point in parentheses, such as 
p{Q), means that it refers to the parameters at the point Q, or is a function of the 
parameters at that point. 

For most cases, however, (standard) notation usually found in the literature 
is used. Vectors are indicated by an arrow A, and unit vectors are indicated by 
a hat a. A dyadic quantity is indicated by a double overline A. Coordinates and 
vectors with a ! (prime) usually refer to the source, and unprimed coordinates will 
usually refer to the receiver or observer. 

1.4.2 Symbols 

These are just a few of the symbols used. Most of them are defined more 
specifically in the text. 
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ttc = cone (half) angle 
(3q = angle of diffraction 
/?s = bistatic angle 

= angle from cone surface 
= kzgz • (f + f') 

£ = permitivity, or a small angle 
£c = f>/a = ellipticity 
r) = y/jije = impedance of material 
6 = theta coordinate, or azimuth angle 
X = wavelength 
fi = permiability 
7T = 3.1415... 

= caustic distance for the diffracted ray 
= caustic distance for the reflected ray 
o = radar cross section 
T = distance along cone surface 
f = tangent vector 
$ = angle from wedge face 
<f>e = <l> for point on edge 
<f>g = (f> for shadow on rim 
= angle from wedge face 
rjje = kcehe • (f + f') 
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u; = 
®1,2 = 


ae = 


a — 
■Aave — 

b = 
6 = 


c 


Ca,Cb,Cc = 

^s,h ~ 

B = 


e = 


= 

e|i = 
= 

E\H^ = 
E\H'^ = 

rC p £r _ 

hi ^Jj — 

jptn Tjm p £r 
£j , ii — 

/ 0 , 1,2 = 
/c0,cl,c2 ~ 

F = 
^ 0 , 1,2 = 


angular frequency 
circular cone fustrum radii 
radius of edge curvature 
radius, semiminor axis 
average spread factor 
semimajor axis 
binormal vector 
speed of light 

partial half-rim corrections 
soft, hard diffraction coefficients 
the dyadic edge diffraction coefficient 
edge vector 

perpendicular direction to plane of incidence 
parallel direction to plane of incidence 
diffracted field 
incident field 

reflected field, radiated field 

field radiated by electric currents 

field radiated by magnetic currents 

equivalent-current half-rim integrals 

equivalent-current corrected half-rim integrals 

transition function 

equivalent-current rim integrals 
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h = cone frustum height 
Ho = zero order Struve function 
Hj = first order Struve function 
_ gjgctric, magnetic currents 
^1,2 = equivalent line current integral 

j = 


Jo = 
Ji = 
k = 
I = 
M = 
n ~ 

A 

n = 
he = 
ho = 
he = 
Qi,2,3,4 = 

Qe — 
Qr = 
r = 

n,2,3,4 = 
re - 
Ts - 


zero order Bessel function 
first order Bessel function 
2?r/ A 
length 

confluent hypergeometric function 

wedge angle parameter 

normal vector 

normal to cone surface 

normal to o-face, endcap 

normal edge curvature 

diffraction points 

point on edge, diffraction point 

reflection point 

distance from coordinate origin 
vector from origin to Qi,2,3,4 
vector from origin to edge 
vector from origin to surface 
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1.4.3 


Tc = 

Ri,2 = 
Ri,2 = 

^s,h ~ 
1 = 
5 = 

S‘ = 

r = 

T = 

'^e,el,e2 ~ 

U = 

U = 

u = 

i; = 

X = 
z = 

Abbreviations 


vector from origin to cone surface 
principal radii of surface curvature 
principal radii of cone surface curvature 
reflection coefficient, soft and hard case 
dyadic reflection coefficient 
distance from diffraction point to receiver 
incident ray direction 
reflected ray direction 

transition function curved-side specular region 
transition function rim-caustic region 
u coordinate in elliptical cylindrical system 
rim-caustic transition function argument 
confluent hypergeometric function 
V coordinate in elliptical cylindrical system 
curved-side specular transition function argument 
z coordinate, or a complex number 


GO = Geometrical Optics 
GTD = Geometrical Theory of Diffraction 
MM = Moment Method 
PO = Physical Optics 
PTD = Physical Theory of Diffraction 
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RB 

RCS 

SB 

UTD 


Reflection Boundary 
Radar Cross Section 
Shadow Boundary 

Uniform Geometrical Theory of Diffraction 
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CHAPTER II 


Theoretical Background 


2.1 Introduction. 

This chapter contains a brief description of the techniques of Geometrical Op- 
tics (GO), Physical Optics (PO), and the Uniform Geometrical Theory of Diffrac- 
tion (UTD). It outlines (without proof) some of the key equations and concepts 
which provide the foundations for this work. More information on these high- 
frequency approximate-solution methods can be found in references [2,3,10,11]. 

In this report, emphasis is placed on the canonical problem of diffraction from 
a curved wedge. For the cone frustum, edge diffraction is a dominant scattering 
mechanism. The canonical problem is important because solutions to more com- 
plex problems can be constructed by the principle of superposition. This will be 
demonstrated in Chapter 3 when the diffraction-point method is discussed. 

2.2 Geometrical Optics (GO) 

The technique of geometrical optics provides a high-frequency approximate 
solution to the incident, reflected, and refracted fields. In our problem, we will be 
only concerned with the reflected field. GO can be obtained from an asymptotic 
solution of Maxwell’s equations, and it corresponds to the leading term of the 
Luneberg-Kline asymptotic expansion for large values of angular frequency [3]. 

According to geometrical optics, electromagnetic waves can be seen as trav- 
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eling along certain lines called rays. The ray path is determined by Fermat’s 
principle. It states that the ray path taken by the light from one point to another 
is an extremum. This is usually a minimum, or the path that requires the least 
amount of time, but not always. From Fermat’s principle, the law of reflection 
and the law of refraction can be obtained. Later we shall see that it can also be 
extended to include the law of edge diffraction. 

In geometrical optics, phase is proportional to the distance along the ray path 
from some reference point, and the amplitude is governed by the conservation of 
power in a tube of rays. Figure 3 shows one such tube of rays. The distances 
Pi and P 2 are from a fixed reference point to the “caustics” , or places where rays 
converge. 

In geometrical optics, however, the fields can not directly be evaluated at the 
caustics. At caustics, the approximation of power conservation in a tube of rays 
is no longer valid. The field near a caustic appears to become infinite as a finite 
amount of power gets squeezed into a vanishing area. The fields near caustics must 
be found by other means. 

2.2.1 GO Reflected Field 


The GO expression for the reflected field (ll| from a smooth conducting convex 
surface with radius of curvature greater than A is (see Figure 4), 


E'{f) = E'(Qr)-R, 




+ »’■) 


-Iks' 


(2.1) 


where 


Qr 

E^iQr) 


S 


r 


reflection point 
incident field at Qr 
distance from Qr to receiver 
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REFERENCE POINT 


; . 


FIELD POINT 



S 


RAY 


Figure 3: An astigmatic tube of rays 
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* UNIT OUTWARD NORMAL VECTOR TO 
THE CONVEX SURFACE AT Op 


A r 

«[r 


A r 

= 6j^X 


Figure 4: Reflection off a convex conducting surface. 
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R 

1 

^a,h 

Pi, 2 


the dyadic reflection coefficient 
Tl 

caustic distances for the reflected rays 


( 2 . 2 ) 

(2.3) 


Pi, 2 


-L + !— 

Pm PgiQr) COS 

2 2 


1± 


Pg(Qr) COs2 e* 

4/2 


J_ _ i. 

A P2 

g cos 2ao 

. Pg{Qr) 


+ 


Pg(Qr)cosg* / J_ _ 

P Vi P\j 


— sin 2oo sin 2uo cos | 


+ 1 


4pg(Qr)cos2 g* 

PR1R2 


(2.4) 


where 


= angle of incidence = cos ^(— s**n.) 


Ri,2 

f 

9 


principal radii of surface curvature at Qr 


1 + 
1 - 


PgiQr) 

PtiQr) 

PgiQr) 

PtiQr) 


COs2 
cos2 6* . 


Note that pg is the radius of surface curvature at Qr in the plane of incidence. 
The plane of incidence contains S*, n, and i. Further, pt is the radius of surface 

A 

curvature at Qr in the plane containing n and 6 where, 


n = normal to the surface at Qr 
i = tangent to the surface at Qr 
h = i X h 
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S* = direction of the incident ray 
S’" = direction of the reflected ray 
= direction of the source from Qr 
r = direction of the receiver from Qr . 

The direction of the reflected ray is defined by the law of reflection which is given 
by 

h-s^ = -h-s* ( 2 . 5 ) 

The point of reflection, Qr, is a point on the surface such that the law of reflection 
is satisfied. At Qr, one finds that 

h’f = h‘f'. ( 2 . 6 ) 

2.3 Physical Optics (PO) 

If the true currents on a scatterer were known, then we would be able to 
find the true scattered field. However, these currents are unknown. The tech- 
nique of physical optics approximates the currents by using the geometrical optics 
currents. For an electrically large conducting scatterer, the surface currents are 
approximately 

_ ^ 2h X H* on the lit surface, and 

J « Jqo = S (2.7) 

0 on the shadowed surface. 

k. 

The PO scattered field is obtained by substituting these approximate currents into 
the radiation integral. 

In this work, PO is used only in an indirect manner. We will use the knowledge 
that PO is known to give good results in the direction of specular scattering. 
References [1,10] provide some discussion on this subject. 
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2.4 Uniform Geometrical Theory of Diffraction (UTD) 

Classical geometrical optics considers scattering due to reflection, but not 
diffraction. Keller [2] extended GO to include diffraction into the approach known 
as the geometrical theory of diffraction (GTD). With the knowledge that wave 
propagation is a local phenomena at high frequencies, it was reasoned that the 
diffracted field, like the reflected field, must also propagate along rays. GTD 
assumes that diffracted rays are produced when incident rays encounter surface 
discontinuites (edges, tips, vertices) or when they graze the surface of a scatterer. 
It is found that various laws of diffraction can result by using these assumptions in 
a generalized form of Fermat’s principle. Experimental results confirm that these 
assumptions are valid for electrically large scatterers. 

In GTD, the initial value of the diffracted rays are obtained by multiplying 
the incident field with a diffraction coefficient at the point of diffraction, analogous 
to the reflection coefficient for a reflected ray. The expression for the diffraction 
coefficient depends on the scattering mechanism (edge diffraction, tip diffraction, 
etc.). The development of new valid diffraction coefficients is presently a major 
area of investigation. The application of GTD is limited only by the availability 
of accurate coefficients. The uniform geometrical theory of diffraction (UTD) is 
an improved version of GTD that uses diffraction coefficients which remain valid 
near the shadow boundaries — where Keller’s original coefficients had failed. 

Two major contributors to the diffracted field considered by UTD are diffrac- 
tion from edges and diffraction from curved surfaces (or creeping wave diffraction). 
In this report, only edge diffraction will be considered. 
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2.4.1 UTD Edge Diffracted Field 


The diffracted field from a semi-infinite perfectly-conducting curved wedge is 
given by the following expression; 

( 2 . 8 ) 

where 


Qe = 

E\Qe) = 


= 


D = 


point of diffraction on edge 
incident field at Qe 
distance from Qe to receiver 
caustic distance for the diffracted ray 
the dyadic edge diffraction coefficient . 


The dyadic edge diffraction coefficient is given by 

3 = -^P'o^oDs - . (2.9) 

The following is the form of the diffraction coefficient which was developed by 
Kouyoumjian and Pathak [3]. The coefficients of a perfectly-conducting curved 


wedge are 

D,h = 7 = 

’ 2nv^2?rA:sin/?o ^ 

+ q: (D’’” + D’’*’)} 

(2.10) 

where 

V 2n 

^ F(ikL*'‘a+($ -$')) 

(2.11) 


=“'‘( 2„ 


(2.12) 


cot 

\ 2n 


(2.13) 


D" = cot f " - 

V 2n 


(2.14) 


21 






The angles /?Qj and are shown in Figure 5, and n is a wedge angle parameter, 
where the wedge angle is (2 — n)7r. The cotangent terms in the diffraction coeffi- 
cient become singular at various reflection and shadow boundaries of the surfaces. 
Table 2.1 and Figure 6 show where each cotangent term becomes singular. The 
function F is a transition function which is defined by 

F{x) = 2jy/xe^^ f dr . (2.15) 

Jy/x 

Figure 7 shows the magnitude and phase of this function. L is a distance parameter, 
where 


/ i.r , \ i,r i.r 

ri,r {Pe + «)Pl Pi .2a 

L ’ = s— ^ ^ sin pq . 


(2.16) 


pr{pY + s){pY + s) 

The parameters pY are the (incident,reflected) radii of curvature in the plane 
containing the (incident, reflected) ray and e; and pj ’2 are the principal radii of 
curvature of the (incident, reflected) wavefront at Qg. The superscripts n,o on 
L denote that the radii of curvature are determined at the respective reflection 
boundaries (2n — 1)7 t — and tt — The reflected field caustic distance in the 
plane of the reflected ray and edge tangent is given by 


where 


1 1 2(n • Ue)(S* • n) 

Pe p\ ae sin^ ^0 


(2.17) 


n = normal to the surface atQ« 
he = normal to edge curvature atQ« 
s* = incident ray direction 
ae = radius of edge curvature 

±/ \ « 2 f 2n7TiV* - x\ 

(x) = 2 cos'^ ( 1 . 


(2.18) 
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Figure 6: Shadow and reflection boundaries for different angles of incidence. 
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Table 2.1: Table showing where cotangent terms are singular. 



The cotangent Is singular 
when 

value of N 
at the boundary 


- 7T, a SB 

surface 4 >e 0 is shadowed 

« 0 

BSBi 

♦* ♦' + 71, a SB 
surface ♦»mr is shadowed 

N‘ » 0 

cot 

= (2n-l)ii-4>', a RB 
reflection from surface 4 >=nTt 

tt 

HiSBI 

4>=7t - a RB 

reflection from surface 41=0 

N" * 0 
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And are integers which most nearly satisfy the equation, 

2nnN^-x=±n (2.19) 

where x = $ ± See Figure 8. 

Law of Edge Diffraction 

The direction of the diffracted rays are defined by the law of edge diffrac- 
tion. This law arises from a generalized statement of Fermat’s principle where the 
diffracted ray path is the minimum distance from source to edge to observation 
point. The law of edge diffraction is given by 

S* • e = S • e (2.20) 

where 


e = unit edge vector 
6* = direction of the incident ray, and 

s = direction of the diffracted ray. 

Note that an incident ray hitting a point on the edge gives rise to a cone of 
diffracted rays (see Figure 5). 

A diffraction point, Qe, is a point on the edge such that the law of diffraction 
is satisfied, or 

- f' • c = f • e (2.21) 

where 


f' = direction of the source from Qe 
r = direction of the receiver from Qe 
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Figure 8: JV"*" and TV as functions of /?q and n. 
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Reasons to Use UTD 


The following are some of the reasons why we choose UTD for our problem: 

1. No exact solution is available for the cone frustum. Exact solutions exist only 
for the most simple shapes (that fit a general coordinate frame). Therefore, 
we have to use an approximate-solution technique. 

2. UTD is accurate for high frequencies. Our problem deals with frequencies 
at the high end of the resonance region and into the optical region. High- 
frequency asymptotic methods, such as UTD, work well at these frequencies. 
Past experiments have shown good agreement with measurements and other 
solutions. The moment method can also provide solutions in the resonance 
region, but computations can become costly with increasing frequency. 

3. It complements low-frequency solutions, such as those obtained by the mo- 
ment method. 

4. It can yield bistatic solutions. 

5. It is easy to apply. The local nature of high-frequency asymptotic approxi- 
mations can resolve a complex scatterer into simpler components. 

6. It is a ray method. Ray methods can provide physical insight into the scatter- 
ing mechanisms. Also due to the continuous nature of physical phenomena, 
the accuracy of a calculated scattering pattern can be readily determined. 
The inaccuracies caused by missing scattering components appear as discon- 
tinuites in the pattern. 

In the next chapter we will show how a solution for the cone frustum can be 
obtained using UTD, and the advantages mentioned above will be demonstrated. 
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CHAPTER III 


The Diffraction-Point Method 


3.1 Introduction. 

The previous chapter briefly discussed the UTD solution to a curved semi- 
infinite wedge. This chapter applies those results to obtain an approximate solution 
to the scattered fields from a cone frustum. The technique used to generate the 
solution is called the “diffraction-point method” . 

3.2 How does the diffraction-point method work? 

We recall that at high frequencies, wave propagation and scattering have a 
very local nature. Ray techniques are therefore applicable. According to the 
geometrical theory of diffraction, diffracted rays occur when incident rays hit geo- 
metric discontinuities on a scatterer, such as the edges. For an edge, the direction 
of these diffracted rays are determined by the law of edge diffraction. 

In order to apply the diffraction-point method, one must first be able to isolate 
the source of these diffracted rays to separate points. The points are sometimes 
called “scattering centers” or “diffraction points” . For the cone frustum, the scat- 
tering centers are located on the edges forming the rim. The location and number 
of points will change depending on the position of source and receiver. Note that 
for the cone frustum, there can be eis many as four diffraction points. 

This local character — our ability to isolate points of diffraction, simplifies 
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the solution to the problem greatly. Since diffraction is a local phenomena, only 
the local geometry in the neighborhood of the scattering center is important in 
determining how the incident energy gets scattered. The diffracted field from a 
scattering center can be approximated with a known solution for a scatterer having 
the same local geometry. The problem of finding the cone frustum’s scattered field 
is now reduced to the problem of finding the diffraction from curved wedges that 
match the geometry at the diffraction points. 

Therefore, the diffraction-point method involves a superposition of semi-infinite 
wedge solutions. The total scattered field is approximated by summing the diffracted 
field contributions from each wedge. This gives a solution that is both fast to com- 
pute and aesthetically pleasing due to its simple form and physical insight of the 
scattering components. 

Again, however, it must be emphasized that the method only works well when 
the local nature of the scattering centers is not violated. The spacing between 
scattering centers should usually be greater than a wavelength. This is true for 
most regions of observation, but in some places the diffraction-point method can 
fail. These difficulties will be discussed later, but happily these problems can be 
corrected using other techniques. 

A short note on the wording: in this report, the “diffraction-point method 
solution” will be defined as the “first-order edge diffraction solution without cor- 
rections”. Where “first-order” means that only the primary interaction of the 
incident ray with the edge is considered. 

Another item worth mentioning is that the previous condition about the 
“wavelength spacing between scattering centers” , should only be taken as a rule-of- 
thumb. The diffraction-point method solution fails gracefully when this condition 
is violated. The solution may be extended to lower frequncies by including higher 
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order terms (multiple ray interactions). 

3.3 History 

The idea of the diffraction-point method is not new. Since exact solutions 
are often difficult to obtain and apply, GTD techniques are usually the most con- 
vienient way to get an approximate high-frequency solution. Some of the people 
who have used this method in the past for the cone frustum include Bectel and 
Ross (4,5,6] and Ryan and Peters [7]. The results of their work show good agree- 
ment with measurement (in the region where the method is valid). Besides the cone 
frustum, this method has been successfully applied to other scatterers, including 
plates and cylinders. 

3.4 Application of the Diffraction-Point Method to the Cone Frustum 

Now let us proceed with how the diffraction-point method can be applied to 
the cone frustum. The procedure consists of the following four steps: 

1. A coordinate system for the scatterer is assigned. 

2. The location of the diffraction points is determined. 

3. The diffracted fields are calculated for each edge corresponding to the local 
conditions at the point of diffraction. This involves finding the incident field 
at the diffraction point, angles of incidence and observation, and the wedge 
angle and curvature. 

4. The diffracted fields from each diffraction point are summed to obtain the 
approximate scattered field. 
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These steps will be shown in greater detail on the following pages. Also given 
are examples of calculated scattering patterns. To simplify things, the examples 
are restricted to a circular cone frustum with source and receiver confined to the 
a plane containing the axis of symmetry. This plane will be called the principal 
plane. (But there is nothing in the theory that limits us to these specific cases.) 

3.4.1 Coordinate System 

For this example, the standard spherical polar coordinate system {r,6,<f>) is 
chosen. When determining the fields (incident or diffracted) for each diffraction 
point, we will need to convert between their own ray-fixed coordinates to this 
system. 

The cone frustum is aligned with its axis along the z-axis and with its center 
at the origin. Figure 9 shows the setup. Here we will define the principal plane 
as being the x-z plane. Later, when considering the orthogonal polarizations of 
incoming and scattered waves, “horizontal” (parallel to the axis) and “vertical” 
(perpendicular to the axis) will be referred to this plane. 

3.4.2 Diffraction Points 

The location of the diffraction points on the frustum’s rims can be found by 
applying the law of edge diffraction. The law of edge diffraction states that 

-f'- e(<^c) = r • e{<f>e) 


or 

(f' + f) . e(<^«) = 0 (3.22) 

where 

f' = X sin &' cos + ii sin 6’ sin 4> z cos (3.23) 

f = I sin 0 cos + y sin z cos 6 (3.24) 
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and for a circular cone frustum, 


e(<^e) = -xsin^e 4- ycos^e . 


(3.25) 


The source and receiver positions are and {d,<f>), respectively. The location 

of the diffraction point on the edge is <j>e. Substituting for r\ f, and c, and solving 
for <f>e yields 


<t>e = tan ^ 


(3.26) 


±(sin $' sin + sin 0 sin <^) 

[ ± (sin 0 ' cos <f)f + sin 0 cos <j>) J 

We find that there are two possible diffraction points associated with each endcap 

rim, or a total of four possible diffraction points for the frustum. In the principal 

plane, the diffraction points are located at = 0 and (j>e = n. These correspond 

to the points Q\, Q21 Qzt and Q4 shown in Figure 10. 

Note, however, that for certain cases all points on the rim satisfy the law of 

edge diffraction (such as for 0 = 0). For these instances, we cannot locate 

isolated diffraction points, and the diffraction-point method cannot be applied. 

These are called “rim caustic regions” and other techniques will have to be used 

to calculate the proper field. More will be said about this later. 


3.4.3 Diffracted Fields from an Edge 


From Chapter 2, the diffracted field from a curved edge corresponding to a 
diffraction point Qe is 

where 

^{Qc) = -Pol3oDs{Qc) - . 

Far Field Approximations 

Several approximations can be made for the far field (see Figure 11). 
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Figure 10: Cone frustum geometry in the principal plane. 
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(3.27) 




5(5 + p ^) 




(3.28) 


and since F{kLu) w 1 as L — ► 00 , 


D 


p J 4 


s.h 


2n\/2nksin Po 


7T 4- ($ — $' 


cot I — ic i) + cot *2n **) + 


2n 


+ R 


5,/l 


(i±L * - t i!i) + cot } 


(3.29) 


The far field expression becomes 

e{ « ^‘(Qe) • \-0okDs{Qc) - ^'^Dh{Qe)\\l • (3.30) 

Incident Field 


Since we assume plane wave illumination, the magnitude of the incident field 
is the same at each edge. However, we still must consider the phase and the 
shadowing of the diffraction points by the structure. The incident field is given by 

E\Qe) = 4)’(Qe)e^’^’‘' (3.31) 


where 


X ! 


Eo {Qe) - 


+ E)^^ if point is illuminated, and 


if in shadow. 


(3.32) 


Depending on the position of the source, different diffraction points are shad- 
owed. In the principal plane, six different regions must be considered. This is 
shown in Figure 12. 
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Figure 12: Regions where different diffraction points are considered. 
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Geometry at Diffraction Points 


The diffraction points for source and receiver in the principal plane are located 
at C?i, Q 2 , Qz, and Q4 as shown in Figure 10. The values of /3q, $, and n 
depend on which diffraction point is considered. The angles $ and are meeisured 
from the o-face. For this problem we define the fiat face (endcap) as the o-face. 
See Figure 13. 

For the principal plane, 


P 


d 


O-e 


ne 


00 


ag sin^ Pq 


• (f + f') 


ai for Qi, Qz 
02 for Q2, <?4 

(3.33) 

X for Qi, Q2 
-X for Qz, Qa 

(3.34) 

for Qi, Q2, Qz, and Q 4 . 

(3.35) 


The value of the wedge angle parameter, n, can be related to the cone angle 
Oc- For a circular cone frustum, 


n- \ ^ ~ ^ 

I I + ^ for Q2, Qi 

where the cone angle is defined as 



for C2 > aj. 


(3.36) 


(3.37) 
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Figure 13: Cone frustum edge geometry. 
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Transformation of Coordinates 


The expression for the diffracted field from an edge uses the coordinates 
(r,/?o,^) which is a ray-fixed coordinate system at the edge. In order to find 
the incident fields and to sum the fields in the next step, it is necessary to convert 
to a common set of coordinates which we defined as (r,^,<^)(see Figure 14). 

Let 


d = r or r' 

ho = unit normal to the endcap (o-face) 
6 = unit binormal vector 
= - (x cos 4>e + y sin 4>e) 

e — unit edge vector, and 

A 

= b X ho * 


In the principal plane, one finds that 


rio — 


b = i 


e = 


+ 2 

for Qu <?3 

(3.38) 

—z 

for Q2> Qa 


+i 

for (? 3 , 

(3.39) 

— X 

for Qi, (?2 


A 

+d> 

for Qi, Qz 

(3.40) 

-4> 

for t? 2 > Qa' 



Since d is the same for both coordinate systems, we can set them equal such that 


A 

d 


X sin 6 cos d> + y sin 6s\n<f>-\- z cos 6 
b sin 00 cos $ + fio sin 0q sin $ + e cos 0q . 


(3.41) 
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Noting that 


A 

ho’ d = sin /?o sin $ 
6 • rf = sin /?o cos $ 


ho • d 

b’d 


tan $ 



and 


00 = cos ^{e-d) 


$ 



00 = — $ X d . 


(3.42) 


(3.43) 


(3.44) 


(3.45) 


The above expressions are best carried out on the computer. Expanding them out 
in terms of 0 and <l> can become a trigonometric mess. But in the principal plane, 


for 4 > — 0 , 


h 

$ 


I for Q 2 > Q 4 
I for Qi, Qz 
+0 for Qi, Q 4 
—0 for Q 21 Q 3 . 


for (j) = n, 


(3.46) 


(3.47) 
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< 


(3.48) 


h = 
4 = 


+<^ 


for Qi, Qz 


-<t> 


for Q2, Qi 


+0 for Q 2 , Q 3 
-0 for Qi,Q 4 - 


3.4.4 The Diffraction-Point Method Solution 


(3.49) 


The diffraction-point method solution is obtained by summing the diffracted 
fields from each edge. The scattered field is approximately 

« Sitd = + £2 + Ei + (3.50) 

where 


Si = E'iQi) ■ 

(3.51) 

Si = S^lQi) ■ 5(Q2)v//(«2)<^’‘’'“ 

(3.52) 

Si = £*(« 3 ) ■ ^(Qz)\J 

(3.53) 

Si = S^lQi) • 5(Q4) . 

(3.54) 


Please keep in mind that Ef = 0 the diffraction point is shadowed or not visible 
(see Figure 12). As an example, the diffraction-point method solution for which 
the source is in Region B and receiver in Region C is 




(3.55) 


which can also be written as 


E^ « <t>Ele 
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(3.56) 


3.4.5 Calculated Principal Plane Patterns Using the Diffraction-Point 
Method 

The patterns shown in Figure 15 and 16 are typical scattering patterns for 
a circular cone frustum in the principal plane which were calculated using the 
diffraction-point method. The horizontal axis is the aspect angle which is measured 
from the 2 -axis. The vertical axis is the radar cross section in decibels above a 
square meter. 

The diffraction-point method solution provides a fairly good approximate so- 
lution, but it doesn’t work everywhere. It fails in the regions near the forward 
scatter direction, the specular direction off the endcaps, and in the specular di- 
rection off the curved side. Since these are directions of significant scattering, we 
cannot ignore these problems. 

The observed failure is due to elements which are missing in the present solu- 
tion. We recall that for some regions, diffracted ray contributions come from the 
entire rim. This is called a “rim caustic” (a system of parallel rays from the rim 
creates a caustic at infinity), and it cannot be handled by the ordinary diffraction- 
point method. From the law of edge diffraction, these caustics occur in the specular 
direction off the endcap (which can involve both front and back rims) and in the 
forward scatter direction. To correct the diffraction-point method solution, the 
contributions from the rims must be considered and added to the solution. In 
Chapter 5, we will show how this can be accomplished. 
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Figure 15: The calculated pattern of a 15*^ cone frustum using the 
diffraction-point method for backscatter. 
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Figure 16: The calculated pattern of a 15® cone frustum using the 
diffraction-point method for bistatic scatter with fixed axial incidence at 0 

degrees. 
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Next we need to explain why the solution fails near the specular direction off 
the curved side. This direction also corresponds to a caustic — since reflected ray 
contributions come from the entire side of the cone frustum. However this caustic 
does not present quite the same problem as for the rim caustic. The diffraction- 
point method solution, to some degree, does take care of some reflection caustics, 
since the diffraction-point method solution does work for a cylinder [l]. Let us take 
a closer look at the diffraction-point method solution near the specular direction 
for the curved side. We note that in Eq. (3.50), as one approaches the specular 
direction, the solution becomes dominated by the diffracted fields coming from the 
two edges forming the reflecting surface. In the region near specular, one finds 
that 

gjfcfi (r+r'l ^ gj*:r 2 (f+r') 


D,^hiQl,2) 00 


and 




+eEl 


DaiQl) Ds{Q2) 


+ 


Dh{Q2 )~ — DhiQi )- — - — 




(3.57) 


The GTD diffraction coefficients for the two dominant edges become infinite in 
the specular direction. “Infinite” because the diffraction coefficients arise from a 
“semi-infinite” wedge solution — the reflected field from a semi-infinite wedge is 
infinite in the far zone. The two dominant edge-diffracted fields have opposite 
signs; each contains the (infinite) diffraction coefficient weighted by the spread 
factor which contains the radius of edge curvature. We note that if the radius of 
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curvature is equal for both edges, as in the case of a cylinder, the total field at 
specular will be correct. The infinities will cancel each other, leaving the correct 
bounded result. However, for the cone frustum, due to the difference in the radius 
of curvature (spread factor) at each edge, this cancellation of infinities does not 
occur, and the solution fails. 

This failure is due to the missing information about the cone’s surface. The 
edge diffraction coefficients have no knowledge of the variation in curvature along 
the cone’s side. This information is necessary in order to determine the correct 
reflected field. If one looks at the situation in terms of rays, there are reflected 
ray contributions to the reflected field from the entire side of the cone (a caustic). 
Each spreads energy differently depending on the radii of curvature at the point of 
reflection. Since this information is not built into the diffraction coefficients, the 
diffracted flteld is not able to go to the reflected field at the specular direction. In 
the next chaf>ter we will look at ways to correct for this failure. 
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CHAPTER IV 


Curved-Surface Specular-Region Corrections 

4.1 Introduction. 

In the preceding chapter, we discussed the failures of the diffraction-point 
method solution. This chapter is devoted to fixing one of them — the failure in the 
specular region for the curved surface. 

As mentioned in Chapter 1, we want to avoid changing back and forth between 
different forms of the solution, such as between GTD and PO. Instead, we would 
like to preserve the GTD form of the solution in all regions. Corrections will be 
accomplished through a modification to the present GTD solution. Two different 
modifications will be presented. One we will call the “spread-factor modification” 
which was shown in a Ph.D. dissertation by T.T. Chu [8] based on similar work by 
Ryan [9]. And the other is a new way that we will call the “equivalent- line-current 
correction” . 

4.2 Physical Optics Result 

In the specular direction, the result based on physical optics gives a good 
approximate value for the scattered field. We will use this result to check our 
solution in later sections. The scattered field (expressed as a radar cross section) 
is given by [4] 

, , A Stt (j3,\ ^ ^ 

a(at specular) = rr-r^ — - cos ( — ) (4.58) 

9A sin^ ttc cos etc V 2 / 
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where /?s is the azimuth bistatic angle (see Figure 10). 
4.3 Spread-Factor Modification 


In Chapter 2, we showed that the diffraction-point solution failed in the spec- 
ular direction due to an incomplete cancellation of infinities. This was caused by 
a difference in the spread factors for the two edges forming the reflecting surface. 
The spread factors contain the caustic distances which differ due to the different 
radii of edge curvature. Chu has shown that by introducing a modified spread 
factor, one can obtain a finite and valid result for the UTD solution in the spec- 
ular direction. The modified spread factor is called the “average” spread factor; 
it is used in place of the actual spread factors for the two edges forming the face 
(n-face). The average spread factor is given by 


lave 



(4.59) 


where 

^ + H-60) 

2 ,4 61) 

3 [ /(Qi) - AQi) J ■ 


Note that the spread factor is no longer local for each edge of the frustum. 

The scattered field with Chu’s modified spread factor can be written in the 
same form as the diffraction-point method solution, Equation (3.50), except that 
the diffraction coefficient is now given by 


^s,h 


-e •'4 


2ny/2irksin0Q 


j^in ^ pro ^ p: 


^ P'avt 

\!aqa 


+ 0") 


(4.62) 
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where 


= cot 1 

1 

+ 

^ 2n 

= cot ( 

1 

1 

^ 2n 

£>’■” = cot 

/tt + ($ + $') 

V 2n 

= cot 1 

^7T - ($ -1- $') 

1 2n 


) 

) 


4.3.1 Calculated Patterns Using the Spread-Factor Modification 

Figure 17 compares the patterns calculated by the original diffraction-point 
method solution with the diffraction-point method solution using the spread-factor 
modification. Note that the spike at the specular direction is now gone. It can be 
verified that the solution using the spread-factor modification gives the physical 
optics result, Eq. (4.58), at the specular direction. However, we can also see that 
the modification has slightly affected the pattern outside of the specular region. 
The modified solution does not exactly go to the diffraction-point method solution 
outside the specular region — it is not a uniform solution. 

The integral in Eq. (4.60) may be viewed as being related to the integral of the 
in-phase currents that give rise to the PO scattered field in the specular direction. 
The lack of a uniform solution is a result of the phcise not being considered in the 
integral for the calculation of the average spread factor. 

Even though the solution is not uniform, the spread-factor modification is 
sufficiently accurate for small cone angles. It has the advantages of being extremely 
easy to apply and fast to compute. But for larger cone angles a better correction is 
needed. Next we will investigate a correction that is based on the use of “equivalent 
line currents”. 
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Figure 17: Bistatic 10° scatter from a 40° cone using (a) the diffraction-point 
method, and (b) the diffraction-point method with spread-factor modification. 
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4.4 GO Equivalent-Line>Current Correction 


The “equivalent current concept” is a technique that can evaluate the fields 
in a caustic region. Recall from Chapter 3 that we have a system of parallel re- 
flected rays off the curved surface of the frustum which creates a caustic at infinity. 
Ordinary GO expressions cannot be used to find this reflected field. However, by 
using “GO equivalent line currents”, the reflected field in the caustic region may 
be determined [1,12]. 

Now one might ask, how can we correct the diffraction-point solution by using 
this technique? This is not an easy question. The approach uses a great deal of 
heuristic argument and physical intuition, but we will show later that it is possible. 
The observations used by K.C. Chiang (ij for the cylinder give us an idea how this 
may be accomplished. 

Chiang had obtained an equivalent-line-current expression for the reflected 
field from a cylinder. She also did an asymptotic evaluation of the physical optics 
solution for an open (no endcaps) cylinder and found that it gave a similar result. 
(And we recall that PO also provides an accurate result in the specular region.) 
Chiang then observed the similarity of the cylinder’s equivalent-line-current solu- 
tion with the portion of the diffraction-point method solution that corresponds to 
the reflection boundary (RB) of the curved surface. And it was then rationalized 
by Chiang that they could be set equal. 

For the cone frustum, we will assume a similar situation where the equivalent- 
line-current solution will yield a PO-type result that is similar in form to the 
RB portion of the diffraction-point method solution (except for the failure in the 
specular region). We cannot justify a direct subsitution for the two solutions, 
but this gives us some insight that it might be possible to provide a “modified” 
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equivalent-line-current solution to replace the RB portion of the diffraction-point 
method solution that fails in the specular region. We will return to the subject of 
correcting the diffraction-point method solution after showing the formulation of 
the equivalent-line-current solution. 

4.4.1 Equivalent-Line-Current Solution Procedure 

The following steps are used to find the equivalent-line-current solution for 
the cone frustum: 

1. First, the GO reflected field over an infinite cylinder is obtained (a two- 
dimensional problem). 

2. Next, the radiated fields from infinite traveling-wave line currents, Z* and 

are determined. 7® is an electric current and 7"* is a magnetic current. 

3. The currents P and are substituited for the cylinder’s surface, and the 
currents 7® and Z”’ are then adjusted so that they give the same field as the 
GO reflected field. In the principal plane, P gives a horizontally polarized 
field, and gives a vertically polarized field. The currents P and 7”* are 
the “GO equivalent line currents” for the cylinder. These are “fictitious” 
currents that vary with the observation position. 

4. The cylinder’s equivalent line currents are broken up into infinitesimal cur- 
rent elements (see Figure 18). The current elements are placed at the equiv- 
alent location of the cone’s reflecting surface, and they are varied according 
to the local radius of curvature for the “equivalent cylinder” at that point. 

5. The radiated fields due to the above currents are evaluated. This involves a 
line integral. The resulting expression is the equivalent-line-current solution. 
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REFLECTED RAYS 



Figure 18: The GO reflected field off the curved side of the cone frustum 
produced by (a) rays, and (b) equivalent currents. 
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It gives the correct value of the reflected fleld off the curved-side of the cone 
frustum. 

The preceding steps will be shown in greater detail on the following pages. 

4.4.2 GO Reflected Field Over an Inflnite Circular Cylinder 

The GO reflected field can be obtained using Equations (2.1) and (2.4) in 
Chapter 2. For a circular cylinder with plane-wave incidence, the expressions 
simplify greatly. For a circular cylinder, R\ = a and R 2 = 00 . The values for pg 
and pt can be determined by using Euler’s equation (see Figure 19 and Figure 20) 
such that 



1 1 2 

— = - cos a 

Pg ® 

(4.63) 


1 1-2 

— = - Sin a . 

Pt a 

(4.64) 

Substituting for Pg^ Pt iii the expressions for p\ 2 yields 


1 

2 { cos^ a + cos^ sin^ a\ 

(4.65) 

1 

1 

1 

a Y cos 6^ j 

1 

P 2 

0. 

(4.66) 

If an angle /? is used, where 

/3 = cos~^ ( 5 ^ • U 2 ) 

(4.67) 

then 

A f A 

^ a 5 ■ n 

2 sin^ /? 

(4.68) 


to 

II 

8 

(4.69) 

Approximations for the far field can be made (see Figure 21) such that 
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Figure 21: Reflected ray off cylinder surface. 
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c> 



( 4 . 71 ) 




gjkre-r ^-jkr' 


And for a two-dimensional cylinder, the spread factor is 


P\P2 


'Pi 


( p 5 + s '-)(/2 + ^’‘) V s ’- • 

Therefore the far-zone GO reflected field is 


(4.72) 


—jkr 




(4.73) 


where the reflection caustic distance p\ is given in (4.68). 

4.4.3 Far-Field Radiation of an Infinite Traveling- Wave Line Source 


Let an infinite traveling-wave electric line current be given by 

_ £^je^-jkucoafi 


(4.74) 


The current flows along the u-axis where l3 is the angle measured between the 
u-axis and the observation direction, s’" (see Figure 22). The far-zone radiated 
electric field |l] is 





(4.75) 


By using duality we can replace 


E -* H 

je ^ jm 




k 


1 

'I 

k 


and we obtain 



TO RECEIVER 
IN FAR ZONE 



Figure 22: Geometry for line source. 
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By use of the plane-wave relation, the fields are 



^ A 

Fk , , 

(4.77) 

if* = 

= (S'x/3)^ 

/JL/V?^ 

1 Sn y/^ 




(4.78) 


= (3 X 3’’)^ 

1 8n yf^ 


If the reference point is moved from the u-axis origin Q to some other point O, 
adjustments will have to be made for the phase. In the far field, one obtains 


5 


r 



in phase, and 
in amplitude 


(4.79) 


where fs is the position vector from the reference point to the u-axis origin. The 
far-zone radiated fields with phase reference at O are as follows: 


j^e 






nm 


^\8n' y/r 


(f X 


Stt 


0 X 

V Stt y r 


a./Aljm yf 


(4.80) 


(4.81) 


(4.82) 


(4.83) 


4.4.4 Cone Frustum Geometry 

To keep the mathematical expressions relatively simple, it is important to 
choose a good coordinate system. For expressions dealing with equivalent currents 
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on the curved surface of an elliptical cone frustum, the elliptical cylindrical coor- 
dinate system {u,v,z) is chosen. Details on elliptical cylindrical coordinates can 
be found in Reference [ 13 ]. 

The cone frustum should be aligned properly to take full advantage of the 
elliptical cylindrical coordinates. The cone axis is the z-axis; the elliptical endcaps 
are surfaces of constant z. The major axis of the ellipse is oriented with the i- 
axis, and the minor axis with the y-axis. The curved surface of the cone then 
corresponds to a surface of constant u. And the equivalent line currents flow on 
the cone surface, along a line of constant v. 

For elliptical cylindrical coordinates, 


that 


where 


I = (7 cosh u) cos V 

(4.84) 

y = (^sinhu) sin v 

(4.85) 

z = 2 . 

(4.86) 

socis, and 6 ( 2 ) be the semiminor axis. 

It can be shown 

a = 7 cosh u 

(4.87) 

b = 7 sinhu 

(4.88) 

a = ai + {z — zi) tan ac 

(4.89) 

b = eca 

(4.90) 

ai = 0 ( 21 ) 

(4.91) 

= 0(22) 

(4.92) 

b 

€c = ellipticity = - 
a 

(4.93) 
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Figure 23: Cone frustum surface geometry 
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(4.98) 


(f X f) X f 
|(f X f) X f 


Ri 2 = principal radii of curvature for the cone. 


R 


C 

1 


a{z) 


tan^ ac + cos v + sin^ 
Ecltan^ Qc (cos^ V + £c sin^ u) + 1] 


(4.99) 


i?2 = oo- (4.100) 

4.4.5 GO Equivalent Line Currents for the Cone Frustum 

The equivalent line currents for a cylinder are obtained by setting the infinite 
cylinder’s GO refiected field equal to the field radiated by the infinite line currents. 
We then solve for the unknown currents. The equivalent line currents for a cylinder 
are 

I "J” 

r = ^.(i»(0).l)y^^^c>*"'"' (4.101) 

= )&-(#‘(0)-l) 

To determine the equivalent line currents for the cone, we replace the position 
vector Ts with fc for the cone frustum, and in pj we replace the radius of curvature 
a with for the cone frustum. may be viewed as being the equivalent cylinder 
curvature at the reflection point located by fc. The equivalent line currents for the 
cone frustum are 

I — ~j~ 

r = ^.(f»(0).l)i/^ — (4.103) 

y k rjy 

= )9.(H‘(0).l)y^r,e-^?^e^*’-"^"' (4.104) 





(4.102) 
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where 



Rlr - he 
2 sin^/? 


4.4.6 Location of Reflection Points on the Cone 


(4.105) 


It is necessary to determine the points of reflection on the surface of the cone 
frustum, so that we will know where to place the equivalent line currents. Let Qr 
be a point of reflection on the surface of the cone. At Qr, the law of reflection is 
satisfied by 

- 3* • n = S'" • n (4.106) 

where n • 3* < 0 (a closed surface), or 

n • (S*^ + 3*) = 0 . (4.107) 


The outward normal to the cone surface is given by 


n = He — 


XEc cos V + y sin v — zec tan Oc 

^ E^ cos^ V + sin^ V + tan** Uc 


(4.108) 


Note that for a given cone, h is only a function of v. In the far zone, the directions 
3^^ and 3* are constant over the surface of the cone. Therefore, for a given 3'" and 
3*, the locus of reflection points is along a line of constant v on the cone’s surface. 
This is a straight line along the side of the cone. The value of v is determined by 
the reflection point condition, such that 


<;x£c cos t; + Cy sin V — <;zSc tan ac —0 


(4.109) 


where 

^ = 5* + S’" = xqx + y<iy + • (4.110) 
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4.4.7 Field Radiated by Equivalent Line Currents 


The equivalent line currents will flow along a line of constant coordinate v 
where the reflection occurs. To And the radiated field, we first use the expression 
for the field radiated by an infinitesimal dipole with uniform current distribution 
Iq and orientation / (see Figure 24) which is given by 


Alts 

The field radiated by a current distribution /(/)/ is found by integration. For a 
straight line source, one obtains 

fh 

E = pjrtksm/3 /(/') dl' (4.112) 

Jil 47T5 


and 

{ r — Tc • r in phase, and 
r in amplitude . 

Therefore the far-zone radiated field is 

- p-jkr fl2 . 

E = PjrjksmlS / . 

47rr i/j 

The far-zone fields due to the equivalent line currents are 


(4.113) 


(4.114) 


Ee 


E., 


m 


^j7?fcsin/?— 

4nr 


0 X r)j k sin /3- 


r^2 

Jti 

r (r')e>*"‘ "dr' 

(4.115) 

—jkr 

4nr 

r /”’(r')e^*"<= "dr' . 
J 

(4.116) 


Now the problem is to evaluate the integral. We might be able to do it 
numerically, but then we do not gain much insight into how the solution behaves; 
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Figure 24: Geometry for radiating current element. 
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besides, it is inefficient in calculation. We would like to get it in closed form — and 
it turns out that we can. 

With a closed form solution, the result is dependent on only the parameters at 
the endpoints of the integration. And we note that the endpoints will correspond 
to the same location as the two diffraction points forming the surface! This gives 
us another clue that it might be possible to combine the equivalent-line-current 
result with the diffraction-point method solution. 

In the expressions for and .E*”, the only parameters that are functions of 
the integration variable r are the radius of curvature and the position vector 
fc (see Figure 25). The integral we want to solve is 


Ju ^ 


(4.117) 


h = 9{v) l’‘ 

Ju 


(4.118) 


where 


?(v) = 


tan^ Qc + cos v + sin^ 


£c [tan^ Oc (cos^ v + sin^ v) + 1] 

By performing a change in the integration variable from t' to z', 


Z - T 


and we obtain 




(4.119) 


(4.120) 


j. = m n 

* 7" 


(4.121) 
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Figure 25: Surface geometry showing f, z, and Fc. 
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where 



d = f + f' = x6x + i/Sy + z6z 

(4.122) 


f z= 6x cos V 4- 6y£c sin v . 

(4.123) 

Now we perform another change of variables from to a, such that 



= a cot Oc + (^l — Ul cot ac) 

(4.124) 

dz' 

— cot occda 

(4.125) 

h 

= cot n y/^eJ^^da 

Z‘f Jai 

(4.126) 

where 

X = k{f + ^iCotttc) • 

(4.127) 

Let 

I2 = / y/ae^^^da . 

(4.128) 


Ol 


The integral I 2 can be evaluated in a couple of different ways. Both give a result 
that uses the confluent hypergeometric function M or U in which case 


J2 — - [ y/ae^^^da + f y/ae^^^da 

Jo Jo 

= -?afe^^“iM(l,2.5,-iXai) 

+ ^o|e^^“2M(l,2.5,-jXa2) . (4.129) 

O 


By using the knowledge that 



(4.130) 
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and adding (4.130) to (4.129) gives 


Jai Ja2 


= afe>^“iU(l,2.5,-iXai) 
-afe^^“2U(l,2.5,-jXo2). 


(4.131) 


Information about the functions M and U is provided in Appendix A. 

The far-zone radiated fields due to the equivalent line currents can now be 
written as 


Ee = 

■ 0\0 ■ (E’(0) • «)| { } 

(4.132) 

Em = 

: ^ xf|§-(£>(0)xf')-S){ } 

(4.133) 

where 



< > " \/r(" • nc)ff(t;)^e^*M^l-‘*i 

(4.134) 


E^-R = e{Rs€^^ + E^RhC^^ 

(4.135) 

{E* X 

f')‘l = EiRhe^j^-E{Rsel. 

(4.136) 

The total field radiated by the equivalent line currents is the sum of fields due to 

both electric and magnetic equivalent currents. It can be separated into parallel 

and perpendicular polarized components such that 



— Ee+ Em 

(4.137) 

E'' 

= (£' ■ ej)e| + (^ . el)el 

(4.138) 
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The following relations can be used in Eq. (4.132) and Eq. (4.133); 


e|i • (^ X f) = ^ • (f X ej:) 


A A T 

r X Cm 


(4.139) 

(4.140) 


We can show that 


• e" - E\Rs { } . 


(4.141) 

(4.142) 


There is no cross polarization. The total field that is radiated by the GO equivalent 
line currents is 


^—Ir . cot ac)!_^^ 


= f;‘(0) • R\^ '—{r . fi,)^(t;)::fi:9£e 

V 47T Z - T 


h . (4.143) 


It can be verified that this expression gives the physical optics result of Equa- 
tion (4.58) in the specular direction. We will call this expression the “equivalent- 
line-current solution” for the cone frustum. 

The equivalent-line-current solution depends only on the conditions at the 
edges of the cone frustum. We can then manipulate the expression into a form 
that resembles a sum of diffracted fields from the two edges. After a great deal of 
work, it can be written in the following form: 


• e\ I 
X, I 




(4.144) 


The points (Qj and Q 2 ) correspond to the same location as the diffraction points 
associated with that surface. The variables (sj and 52 ) the distances from the 
points {Q\ and Q 2 ) to the observation point. The expression contains terms where 
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the incident field at the edge is multiplied by a factor called the “physical-optics- 
type” diffraction coefficient for that edge and a local reflection spread factor. We 
call it a “physical-optics-type” diffraction coefficient because its form is similar to 


the physical optics diffraction coefficient for a half-plane. The PO-type diffraction 


coefficients are given by 
,PO- 


D 

D 


3jh 

PO^ 

5,/l 


= -R 


-7^ 

e ^ 4 


sin/? 


-^2irk cos /? + cos /?' 


- 1 - 
e ^4 




where 


X = kcoi ac{f • (f 4- f')) 


(4.145) 

(4.146) 


(4.147) 


and the function T is called a “transition function” . 


Transition Functions 


The function T, that we isolated from the equivalent-line- current solution, 
has different forms depending on how the integration in Equation (4.128) was 
performed. If it was evaluated as in Equation (4.129), we obtain the transition 
function, 

T(0) = Ti(z) = ^2M(1,2.5 ,-z) . (4.148) 

Determining the integral bls in Equation (4.131) results in 

T(z) = T 2 (z) = -zU(l,2.5,-z) (4.149) 

where 

U(l,2.5, -z) = -|m(1,2.S, - 2) + . (4.150) 

Physically, T^ may be viewed as a result due to a finite cone frustum, and T2 due 
to a semi-infinite cone frustum. Mathematically, both give the same result when 
both edges of the reflecting surface are considered. 
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The confluent hypergeometric functions, M and U, have series and large ar- 
gument forms which can be calculated rather quickly on a digital computer (see 
Appendix A). By using the first few terms of the series and large argument forms, 
the small and large argument forms of the transition function are the following: 

for l^l 1, 


Ti(z) « 


/2z 2z^ z^\ 

(4.151) 

T2(^) « 

e~^ 

f y/ir 2z 2z^ z^\ 

V2V^^ 3 5 7 J 

(4.152) 

for \z\ » 1, 




TiW ~ 

1 + 

_ i L 

2 yfz 2z 4 z 2 

(4.153) 

T2(j) ~ 

1 - 

1 1 
2z 42^ 

(4.154) 


The magnitude and phase of Ti and T 2 are shown in Figures 26 and 27 for values 
of the argument z = jx. The transition functions have their greatest effect in the 
specular region (where the value of the argument is small). The transition functions 
have little effect outside of the specular region. Both Ti and T 2 approach unity 
for large values of the argument. 

4.4.8 Modified Equivalent-Line- Current Solution 


Now we will return to the question of how to correct the diffraction-point 
method solution in the specular region. In the equivalent-line- current solution, 
we obtained a physical optics type diffraction coefficient which is given by 

i!!LL_^T(±iXa,) . 

’ \/ 27 rA: cos/? -f- cos/?' 
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It resembles the PO diffraction coefficient of a half-plane multiplied by the function 
T. We mentioned earlier that Chiang's asymptotic evaluation of the PO solution for 
the open cylinder yielded the cylinder’s equivalent-line-current solution. Therefore, 
for the cone frustum, we shall assume that cone frustum’s equivalent-line-current 
solution is the asymptotic PO solution of an open cone frustum. 

Let us now proceed with the following line of reasoning towards a “modified” 
equivalent-line-current solution that may be used to correct the diffraction-point 
method solution. If the exact currents on the cone frustum were known, we could 
determine the exact scattered fields. But we do not know the exact currents, so 
we must rely on the approximations provided by the methods of physical optics 
and GTD. We know that physical optics provides an accurate representation of 
the scattered fields in the specular region, but it does not work well outside the 
specular region due to errors in the edge currents. GTD, on the other hand, 
provides accurate edge currents, but it fails in the specular region (for the cone 
frustum) due to the lack of information about the reflecting surface. Somehow 
there must be a way to combine the desirable characteristics of both PO and GTD 
to provide an overall accurate representation of the scattered field. This is similar 
to the philosophy of the physical theory of diffraction (PTD) where a correction, 
called the fringe current, is added to the PO solution. 

Let us further investigate the PO-type diffraction coefficient in the equivalent- 
line-current solution. We note that the information on the surface curvature is 
contained within the transition function T. The transition function is active only 
in the specular region. Outside the specular region, we have the PO diffraction 
coefficient of a half-plane — which is known to be inacurrate. 

Let us confine our observation to the plane perpendicular to the edge (see 
Figure 28). The angle 0 is converted to ip which is measured from the curved 
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Figure 28; Edge geometry showing the angles 0, and 



surface. Outside the specular region, one finds that 


T (outside specular region) —*■ 1 
’ v^27rik cos if + cos 


Let’s look at the following factor: 




— sinV' 


(4.155) 

(4.156) 

(4.157) 


cos xjf + cos i>^ ' 

By using the appropriate trigonometric identities, this factor can be rewritten in 
the form 


. (4.158) 

This expression contains terms which indicate a sum of “difference currents” on a 
flat plate, which is what one obtains using physical optics. Note the negative signs 
associated with two of the cotangent terms. If the fringe current (the correction) 
is added to the PO result, we obtain the GTD result. This can be accomplished 
by replacing the negative cotangent terms with plus signs. Therefore the modified 
expression is 




+ cot 4- cot I . (4.159) 

If the edge (wedge angle parameter n = 2) is opened up to form a wedge, the 
cotangent terms are again modified such that 

— I-") 
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Note that this does not affect the result in the specular region, since it is relatively 
insensitive to the wedge angle there. If only the outer curved surface information 
is retained, one obtains 

'‘■•■I 


The dominant term is the reflection boundary term, so 

-{xP + tP> 




(4.162) 


The PO-type diffraction coefficient after these modifications is 

nroT = ji ,= 0 ±cot ( lSlSl±m T{±JXa,). 


(4.163) 


Converting the angle V’, which is measured from the curved surface, to the angle 
which is measured from the endcap face, gives 

+ ($ + $' 


— C ^4 1 

= Rsh ‘ cot 


s,h 


’’ y/2nk 2n 




2n 


T(±jXa*). 


(4.164) 


The expression is made three-dimensional again by inserting the factor sin/0o- The 
modified PO-type diffraction coefficient is 

= R, ^ . cot ( ?L±iJ±iO'\ r(±jXa,). (4.165) 

''^2n\^sml3o \ 2n y v ./ «; v / 

Now we will make another change to the equivalent-line-current solution. The 
reflection caustic distance, is only valid at the specular angle. We will replace 
this with the caustic distance for the diffracted ray, which will give the correct 
spreading of the rays. (At the specular angle, pj = p*^.) 

To summarize, we make the following substitutions in the equivalent-line- 
current expression: 


sin (3 


-1 


cos 13 -t- cos /?' 


2n sin 


cot 
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and 


Pi 



(4.167) 


This converts the PO-type result to a UTD-type result. The modified equivalent- 
line-current expression is 




«1 

,-jks2 

52 


where 


= Rs,h 


-c 


= R, 


2n\/27rfcsin/?o 

V*' 

-C 


cot 






and 


(4.168) 


(4.169) 

(4.170) 


X = A:cotac[tan^ac(cos^v -H e^sinw) -I- 1]5 

• .o + / 2n7T - ($ -I- $') \ 1 

•sm/?ocos( — - — j|2cos( ^ . (4.171) 

Note that Equation (4.169) without the transition function is the far-field n-face 
RB term of the edge diffraction coefficient given in Equation (2.13). 

Figure 29 shows the curved surface RB portion of the diffraction-point method 
solution, and the equivalent-line-current solution. The equivalent-line-current so- 
lution provides a good result in the specular region which includes the main beam 
and first sidelobe. Figure 30 compares the field radiated by the equivalent line 
currents with the field radiated by the modified equivalent line currents. Note the 
good agreement in the main beam and first sidelobe. 

The modified equivalent-line-current result using different transition functions 
is shown in Figure 31. Note that for observation points above the surface, where 
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both edges are present, the pattern is identical. However, below the surface (angles 
greater than 130 degrees), where only one edge contributes to the field, the pattern 
is different. This difference is due to the different natures of the transition functions 
used (see Figures 26 and 27). The transition function Ti oscillates and dampens 
slowly to 1 outside the specular region, while the function T 2 goes to 1 rather 
quickly. In the modified equivalent-line-current solution, when both points on the 
edge contribute to the field, the same result is obtained, but when only one point 
is present the differences become apparent. 

Which transition function should we use? For a finite cone frustum, we can 
use either Tj or T 2 . The small differences in the scattered field below the surface 
can be neglected, except at low frequencies where T 2 should be used (since the 
value of the function argument remains small). The transition region of T 2 is much 
narrower than Ti. Also for a semi-infinite cone frustum, T 2 should be used, since 
without the other edge, we would see a oscillatory nature if Ti was used. Overall, 
theoretically' speaking, ^F '2 produces a better result it provides a more XJTO— 
type” solution. However, from a practical standpoint, Tj is better for calculations 
on the computer, since we do not have to deal with an infinity when the argument 
becomes zero. 

4.4.9 The Diffraction-Point Method Solution with Equivalent-Line- 
Current Correction 


By replacing the curved-surface reflection boundary term of the original GTD 
diffraction coefficient with Equation (4.169) from the modified equivalent-line- 
current solution, we obtain a new diffraction coefficient. The new diffraction coef- 
ficient is 

S’" = -3o3oO? - i'icj 
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where 


= 


-e^4 J f n + — $') \ /tt — $ 


2nv/27rA:sin/?o 


l»'(^) 




+ ^S,/l 


cot 


7T + ($ + $') 


2n 


T(±jXae) + cot 




2n 

- ($ + ^0 

2n 


or 


where 



(4.172) 

Ra,h = Tl . 

(4.173) 

;ion is either 


T(«) = Ti(z) = ^2M(1,2.5,-«) 

(4.174) 

T(z) = T2{z) = -2U(1,2.5,-2) 

(4.175) 

T (outside specular region) —* 1 . 

(4.176) 


The diffraction-point method solution with the equivalent-line-current correction 
for the curved-surface specular region is 

S’ ^ = Ef + + St + Bi (4.177) 


where 


— yifcr 

Sf = .g'(Q,) • 3*(Q,)sJpHQi)e>l‘’-'i ’‘— 

St = £‘((32) ■ s 


—jkr 


St = S'm ■ 

Si = S'{Q,) ^ — 


(4.178) 

(4.179) 

(4.180) 

(4.181) 
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Note how the original diffraction-point method form of the expression is preserved 
by the use of the transition function in the diffraction coefficient. 

Figures 32, 33, 34, and 35 demonstrate for various cone angles how the tran- 
sition function can correct the original diffraction-point method solution in the 
specular region. Note that it is a uniform solution. The diffraction-point method 
solution remains unaffected outside the specular region. Additional information 
on the calculated patterns is contained in the following section. One item that is 
worth noting is the relative CPU times required to generate the patterns. The cor- 
rected diffraction-point method solution takes longer to calculate, but it is roughly 
of the same order of magnitude. For Figure 35, 0.04 minutes versus 0.07 minutes 
on a VAX 11/780 computer with a floating point accelerator. 

4.5 Background Information on Calculated Patterns. 

This section contains information on the patterns for the scattered field calcu- 
lated in this chapter. The patterns were generated using a VAX 11/780 computer. 
All patterns are for a circular cone frustum in the principal plane for the horizontal 
(hard) polarization case only (see Section 3.4.1). Only one polarization is studied, 
since the plots in this section are for demonstration purposes only. Actual com- 
parisons with measurments and the moment method calculations will be shown in 
Chapters 6 and 7. 

The vertical axis of the plot is the radar cross section in decibels above a square 
meter. The horizonal axis is the angle d (azimuth) for the receiver position in the 
principal plane (see Figure 10). Points were calculated every one degree. The 
following table lists information about the cone frustum dimensions, frequency, 
bistatic angle, and CPU times for each plot. The symbols used in the Table 4.1 
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40 


ANGLE 


0.0 OEGREES) 



Fieure 32: Backscatter from a 15° cone frustum using (a) the diffraction-point 
method, and (b) the diffraction-point method with equivalent-line-current 

correction. 



(b) 


ANGLE ( 0.0 DEGREES) 


Figure 33; Backscatter from a 30° cone frustum using (a) the diffraction-point 
method, and (b) the diffraction-point method with equivalent-line-current 

correction. 
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21 


61 
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80 


RNGLE ( 10 


.0 DEGREES) 



iiure 35; Bistatic 10° scatter from a 40° cone frostura osinS (») .‘k' 
Ktion-point method, and (b) the dilTraction-point method with 
equivalent-line-current correction. 


refer to the following: 


/ = frequency (GHz) 
h = cone frustum height (inches) 

) ®2 = cone frustum radii (inches) 

ttc = cone angle (degress) 

/?a = bistatic angle (degrees) 

Tcpu = CPU time required (minutes) 
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Table 4.1: Table of dimensions, angles, and CPU times for calculated patterns 

shown in Chapter 4. 


Figure 

n 

D 


02 

ote 

m 

Tepti 

17(a) 

6.05 

3.75 

1.00 

4.15 

40 

10 

0.04 

17(b) 

6.05 

3.75 

1.00 

4.15 

40 

10 

0.05 

29(a) 

6.05 

3.75 

1.00 

2.00 

15 

0 

0.04 

29(b) 

6.05 

3.75 

1.00 

2.00 

15 

0 

0.07 

30(a) 

6.05 

3.75 

1.00 

2.00 

15 

0 

0.07 

30(b) 

6.05 

3.75 

1.00 

2.00 

15 

0 

0.07 

31(a) 

6.05 

3.75 

1.00 

4.15 

40 

10 

0.06 

31(b) 

6.05 

3.75 

1.00 

4.15 

40 

10 

0.06 

32(a) 

6.05 

3.75 

1.00 

2.00 

15 

0 

0.04 

32(b) 

6.05 

3.75 

1.00 

2.00 

15 

□ 

0.07 

33(a) 

6.05 

3.75 

1.00 

3.17 

30 

0 

0.04 

33(b) 

6.05 

3.75 

1.00 

3.17 

30 

0 

0.08 

34(a) 

6.05 

3.75 

1.00 

4.75 

45 

0 

0.04 

34(b) 

6.05 

3.75 

1.00 

4.75 

45 

0 

0.07 

35(a) 

6.05 

3.75 

1.00 

4.15 

40 

10 

0.04 

35(b) 

6.05 

3.75 

1.00 

4.15 

40 

10 

0.07 
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CHAPTER V 


Rim Caustic Corrections 


5.1 Introduction. 

In the last chapter, we developed a uniform correction to the diffraction-point 
method solution for the frustum’s curved-surface specular region. One can easily 
appreciate the use of a transition function to preserve the original form of the 
diffraction-point method solution. Numerical integrations are avoided, and the 
transition functions have large argument forms which greatly reduce calculation 
times. 

In this chapter we will continue in the same tradition — to develop a uniform 
correction using transition functions for the rim-caustic regions. Fortunately, a 
great deal of the work on this problem hzis already been done. Ryan and Peters [7] 
had developed a solution in the rim-caustic regions based on the fields radiated by 
equivalent edge currents on the rim. Since this result goes to the diffraction-point 
solution outside the caustic region, it is a uniform solution. Recently Chiang [l] 
extended the work of Ryan and Peters by using equivalent currents that use the 
"stripping concept” [15,8] in addition to using small-argument approximations in 
the radiation integral. Chiang’s solution provides an accurate result in the rim- 
caustic regions that are associated with a specular or forward-scatter direction of 
the endcaps. This solution is called the “Bessel-Struve function extension”. It, 
also, is a uniform solution. 
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Uniform solutions for rim caustics are available. However, in the literature, 
the solutions are not often expressed in the transition function form that we desire. 
Also the partial illumination and partial observation of a rim, as in the case of the 
cone frustum, is not considered. These details will be investigated in this chapter. 
First we will look at the equivalent-edge-current solution of Ryan and Peters. 

5.2 Ryan/Peters Equivalent-Edge-Current Solution 


The formulation of the equivalent edge currents is similar to the concept used 
in the last chapter. It will not be discussed here, since much is available in the 
present literature. Some suggested references are [14] and [7j. The equivalent edge 
currents used by Ryan and Peters are 


where 


and 




(5.182) 

(5.183) 

(5.184) 

(5.185) 

(5.186) 

(5.187) 


The equivalent electric current is I‘ and the equivalent magnetic current is 
(see Figure 36). These currents are distributed on the illuminated edges of the 
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Figure 36: Rim geometry. 
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cone frustum and the radiated fields are evaluated. For simplicity, we will restict 
ourselves to a circular cone frustum (of radius a) and observation in the principal 
plane. The far>zone radiated fields due to the electric and magnetic equivalent- 
edge-current loops ue 

^ (5.188) 


27T 


ET = 


El, = 


rptn _ 
- 


ar/e-J^ f 

2ir Qtn 

2tt r Jq 

1 siny^o* 


2)T Qt 

2n r Jq 

1 sin^So 

at] cos B t 

-jkr r2n 


(e • ^*(0)) cos(^ - (5.190) 


2ir 


(5.191) 


The diffracted field is the sum of the fields radiated from both electric and magnetic 
equivalent currents, such that 


>d 

rim 

BEg + ^E^ 

(5.192) 

= 

E% + E^ 

(5.193) 

Ee = 

El + E^. 

(5.194) 


To simplify the the integration in (5.188)-(5.191), Ryan and Peters used the 
approximation that the diffraction coefficients remained relatively constant with 
respect to the variable of integration, so that 


G«.m 

— — « constant in 
8in/?o 


(5.195) 


Let 


= kahe(<^e) ' (f + r') 


(5.196) 
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= jka[sin 9* cos(^^ — + sin cos(^ — <l>e)] 

Tfe = kzez-{r + f'). 

Using the previously mentioned approximation and the fact that 

f sin (f>e cos d<f>e = 0 

Jo 

the expressions for the radiated fields become 


(5.197) 

(5.198) 


(5.199) 




a e~jkr 


i — E\ I -G^ f COS^ <!>e^'^‘d<t>e 

* 27rsin/?o r * \ Jo 

+G^ cos e cos sin^ <t>eeJ^*d<f>c | (5.200) 


Eg = —— — Eg | g * cos ^ cos ^^ f Bin^ ^ee^^*d<i>e 

^ 2jrsin/3o ^ ^ I Jo 


+G 


J cos^ 


where 


rl>e = U cos <f>e . 


(5.201) 


(5.202) 


Let the following integrals be defined by: 


FdU) = 

f cos^ d<j>e 

Jo 

(5.203) 

hW) = 

t — 

J * cos^ <j>ee^^^d<f>e 

(5.204) 

nw) = 

3ff 

1 * cos^ (j>ee^^‘d<f>e 

•'f 

(5.205) 

F2(U) = 

cos ^ cos f siv^ <l>ce^^^d<l>e 

Jo 

(5.206) 
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( 5 . 207 ) 


hW) 

fm 



cob9cos0 





These integrals can be evaluated in closed form as follows: 


h(v) = »[(jo(t7)-^)+>(Ho(tr)-5^)j 
n(V) = >r[(jo(l')-^)-y(Ho(t')-5^)] 

= co,5co,5',(iM^,H^) 
rm = e<.5cc*v(iM.,H^). 

The expressions for the radiated fields are 


The above can be rearranged into 


( 5 . 208 ) 

( 5 . 209 ) 

( 5 . 210 ) 

( 5 . 211 ) 

( 5 . 212 ) 

( 5 . 213 ) 

( 5 . 214 ) 


103 


= -Mr 

(5.215) 

r B5{^(B. + F,) + B(F,-B,)}. 

(5.216) 

Breaking the integrals apart into 


Fi)!!) = /i(ti) + /r(ti) 

(5.217) 

FiiU) = h{V) + n{V) 

(5.218) 

and further manipulation gives 



+^(/r - /J) - B(ll + /J)} (5-219) 

+A(f\ + H) + B(/; - /j’)} . (5.220) 


Now let us associate a half-rim integral with the points at = 0, tt which are the 
stationary points for large V’e^ (These are the diffraction points on the rim for the 
principal plane. See Section 3.4,2.) In the above equation, the field contributions 
from the half-rim associated with = 0 are 


E\4> 


27t sin/3o 


g-j*r 


r 




(5.221) 


104 


(5.222) 


Sf = . (5.223) 

Comparing this with the expression for the edge-diffracted field from a point on 
the rim, 

Bi = E^(Q.) ■ (5.224) 

we obtain the following new diffraction coefficients: 

^[Qt) = -Pq^qDs{Qb) - ^‘^Dh(Qe) 

^ \ (^T.i(Cf)-BT,2(I/)) (5.225) 

2nv27rfcsm/?o 

^ + BT,i{U)) (5.226) 

2nv27rfcsinpo 

where 

T.,(l/) = l(MU) - (5.227) 

T,2(l') = l(hW) + (5.228) 

and 

f^(Qe) = A:o*ne • (f + r ') . (5.229) 

The functions Tgi and Tj 2 are the transition functions through the rim caustic 
region based on the Ryan/Peters equivalent-edge-current solution. The transition 
functions have small and large argument forms which were obtained using the 
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series and asymptotic forms of the Bessel and Struve functions, such that 


for |t/| < 1 

T.i(y) 

for |t/| > 1 

T,2(V) 




• {H* ■ ^ ^ t)] 

(‘ - ^)1 } 




nU 


fir 3t/2 ^ 


.21/ [ ( 2U^\ ,( 




39 


J45 


8f/ 128C/2 i024f/3 

3 Jl5 \ 


— cos Q cos Q' 


\ u ' su^ USU^J 




(5.232) 


l-^ + 


39 


J45 


SU 128C/2 i024t/5 

3 Jl5 \ 


+ cos 6 cos 0' 


' (-^ 

V u ' 8 U^ 12SU^J 

[cosOcosliy (l + i) 


2 

1/2 


(5.233) 
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The total field radiated by the equivalent edge currents on the rim may be 
written in the difTraction-point form. 

= ^ + ^ (5-234) 

where 

(5.235) 

Et = -E’(<?3) • B(<?3) ■ (5.236) 

Note that Q\ and are the two diffraction points on the rim. 

5.3 Bessel-Struve Function Extension 

One problem associated with the Ryan/Peters result is that the approximation 
(5.195) becomes invalid for a rim caustic associated with an endcap’s specular or 
forward-scatter direction. The diffraction coefficients can become infinite in the 
specular or forward-scatter direction. 

It was also discovered that by evaluating the integrals in (5.188)-(5.191) nu- 
merically, the result is twice the expected result in the specular rim-caustic regions. 
The explaination for this is that the reflecting surface is evaluated twice when the 
integral is performed around the rim from 0 to 27t. To overcome this difficulty, a 
modified edge vector e* is used. It is defined by 

e* = p{e • p) (5.237) 


where 



no = normal to the endcap 


(5.239) 

f' = direction of source, and (5.240) 

f = direction of receiver. (5.241) 

By using the modified edge vector, only the perpendicular components (to the 
plane of incidence) of the equivalent edge currents are considered. This is called 
the “stripping concept” . The equivalent edge currents with this modification are 
given by 


II = 

^ r;Asin|Sosin/?Q 




(5.242) 




(5.243) 


The modified equivalent-edge-current solution is obtain in the same way as 
before. To simplify the evaluation of the integral, Chiang used the following small- 
argument approximations: 


6 « $' — € ,in the specular region, or 
0 « 7T — , in the forward scatter region. 


This resulted in 


and the integral 


G* 


m 


sin^osin/?Q ^ ccos^e 


P'oiU) = 

Jo 


(5.244) 

(5.245) 


(5.246) 


(5.247) 
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which can be evaluated by using 


/o(I7) = f 

•'-f 

(5.248) 

= <r J-H,(t;)+M(l')] 

(5.249) 

.to 


f^iU) = / cos 

•'f 

(5.250) 

= X -- + Hi{U)+jJi{U) . 

7T 

(5.251) 


The radiated fields using this approximation, give the physical optics (projected 
area) result in both the specular and forward scatter directions. Associating a half- 


rim with a diffraction point, the expression of the edge diffracted field becomes 

Ej = & . 

2 X r ' 


Chaing’s solution is a uniform solution. It can be written in the diffraction-point 
form by using a transition function in the diffraction coefficient, such that 


EL = Bi + Ei 


(5.253) 


where 

E{ = eHQi) (5.2M) 

H = E'(Qi) ■ . (5.255) 

Again, Qi and are the two diffraction points on the rim. The diffraction 
coefficient is 

^{Qc) = -0O^oDs{Qc) - ^'^£>fc(Qe) 
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(5.256) 


where 


^s,h ~ 


-e->f 


2 nv/ 27 rA: sin /?o 



(5.257) 


U{Qe) = kaehe-{f + f') 


(5.258) 


The transition function is Te(l^). The following are the small and large argument 


forms: 


for \U\ < 1 



(5.259) 


for \U\ » 1 


T«(t/) 




15 

128t/2 


jT05 

1024C/3 


(5.260) 


The magnitude and phase of the the transition function T; as a function of the 
argument U is shown in Figure 37. 

5.4 Rim Shadowing 

The previous equivalent edge current solutions assumed no shadowing or 
blockage of the rim from observation by the structure. In the actual situation, 
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however, we must consider partial rims — rims which are partially illuminated or 
visible. In order to know the size of these partial rims, we must first determine 
where the shadows are located. 

The expressions for the surface normal and the direction of the source are 
known. For a circular cone frustum with C 2 > oi, these are given by 

he = X cos ac cos + y cos Oc sin ^3 + 5 sin Oc (5.261) 

f' = X sin o' cos <i>' + y sin O' sin <f>' z cos O' . (5.262) 

At the boundary between the illuminated and shadowed portions of the cone frus- 
tum (the terminator), one obtains 


or 


where 


hc-r' = 0 

A cos <^3 + J9 sin ^3 -h C = 0 

A = cos ac sin O' cos 4>' 
B = cos ttc sin o' sin <f>' 

C = sin Qc cos o ' . 


(5.263) 

(5.264) 

(5.265) 

(5.266) 

(5.267) 


The unknown that we want to find is <f>a. Let us make the following substitution 
for cos ^3 and sin ^ 3 : 


cos <i>a = 



(5.268) 

gj4>» 

(5.269) 

2 


■ X 
2 

(5.270) 
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sin <!>i 


fj4>$ _ 


(5.271) 


2i 



Equation (5.264) becomes 


(5.272) 


[B + jA)x^ + j2Cx + (-B + }A) = 0 . (5.273) 

We can then solve for x by using the quadratic formula. Solving for x yields 


Re{x} — 


sin Q!c cos ij>' cos 0^ T sin 4>*\/ cos^ Oc — cos^ 0’ 
cos ae sin 0* 


(5.274) 


7m{x} = 


sin Oc sin <f>' cos 0* ^ cos cos^ olq — cos^ 0' 
cos ttc sin 0' 


(5.275) 


The values of <f>a that correspond to the boundary between shadowed and illumi- 
nated portions on the cone frustum are 


^31,32 = tan 


-1 



(5.276) 


= tan 


_l f sin Oc sin cos 0' ^ cos 4>’y/coB^ ae — cos^ 0' 
sin Oc cos <t>* cos 0' ^ sin (^y/ cos^ ae — cos^ 0' 


(5.277) 


The above equation applies only for the case where cos^ ac > cos^ 0'', otherwise, 
the rim would be completely illuminated or all in the shadow. 

In order to know how the shadowed rim appears, we need to know the region 
where the source is located (see Figures 12 and 38). The results of this section 
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apply equally well to finding the portion of the rim that is blocked from view by 
the cone structure. 

5.5 Partial Rim Corrections 

In the equivalent-edge-current solutions, we associated a half-rim with each 
edge diffraction point. As a result, there are half-rim integrals in the transition 
functions of the diffraction coefficients. If the shadowing of the rim is considered, 
we will have partial half-rim integrals. In this section we will provide corrections 
to the half-rim integrals by subtracting the portion that is shadowed or not visible. 
Let 


MU) 

= J 

P cos 

(5.278) 

MV) 

= J 

f ^ cos^ <l>ee^^^d<f>e — Cg{U) 

“f 

(5.279) 

MV) 

n 

= cos 6 cos 6' j s\l? — Cc[U) . 

•'-f 

(5.280) 

The correction terms 

are given by C^, C^, and Cq (see Figure 39). 

If the center 


of the half-rim is missing the correction terms are 


Ca{U) = / cos<f>ee^^*d4>e 
hi 

(5.281) 

Cb{U) = f*\os'^4>ce^^*d<f>e 
hi 

(5.282) 

Cc{U) = cos ^ cos f SIX? d>e^'^‘d<i>e 

J 

(5.283) 
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0,4 ^-0 
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\ 



a. 

^“T 

PARTIAL HALF- RIM 
^^ENDS MISSING 
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/ 


\ 


\ 





j. -JE- 

9 * 2 


Figure 39: Two types of partial half-rims. 


116 



where — | < 0, 0 < ^2 < f » and |^i| = \<t> 2 \. If the ends of the half-rim are 

missing the correction terms are 


Ca{U) 

csm 

Cciu) 


2 / cos 4>e^^*d<l>e 

■'-f 

2 / cos* 

•'-f 

2 cos $co8^ j ^ sin* <l>ee^^*d<f>e . 


(5.284) 

(5.285) 

(5.286) 


These integrals may be evaluated for small values of \U\ by the following series 
form obtained by expanding the exponential. It converges rapidly due to the it! in 
the denominator of each term. This is sufficient, since for |l/| greater than about 
8, we are out of the transition region, and the correction would not be necessary. 

rrik t 

COS(^eC^^*d(^e = / cos^+'^ed^e (5.287) 

ifc=0 

Tj^ r 

cos^ <f>ee^'^*d<t>c = X^y*^-7T / cos*'*'* ^ed^e (5.288) 

*=0 ■ •' 

TJ^ t 

sin*^ee^'^*d(^e = / cos* ^ed^e - ^^(Cf) (5.289) 

*=0 

and 

J cos” 4>ed<t>e = ^ cos”"^ 4>e sin 4>e + j COs”“* <t>ed(f>e . (5.290) 

The number of terms required in the series increases as |17| increases. 

Tables 5.1 and 5.2 show the diffraction points associated with a partial half- 
rim, the type of partial half-rim, and operations to determine the limits of inte- 
gration <i>i^2 for various regions of observation and illumination. The symbols used 
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in the table refer to the following: 


ii* = illumination region 
= observation region 




Lq2 = 


Lq3 = 


I'QA = 


L 


cen — 


^end = 


> = 


< = 



Q\ considered 
Ql not considered 

Q2 considered 
Q2 not considered 

considered 
Q3 not considered 

Q4 considered 
Q\ not considered 

partial half-rim with center missing 
complete half-rim 

partial half-rim with ends missing 
complete half-rim 

4 >al,s 2 with greater |<^,2 ~ are ^1,2 


<^ai,s2 with smaller \ 4 >s 2 ~ are <^1^2 
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A* 


A*" 


^ 1 ^8l,s2 of source region are 

0 — 

k. 

{ 1 ^«l,«2 of observation region are ^ 1^2 

0 - 


The following is a step-by-step procedure to find the correction terms: 

1. Determine the regions of the source and receiver (see Figure 12). 

2. Consult Tables 5.1 and 5.2. Check if a partial half-rim is present, and the 
diffraction point associated with each partial half-rim. If no partial half-rim 
is present or if |l/| > 8, then no correction is necessary. 

3. If a partial half-rim is present, find the values of <l>gi^g2 given by Equa- 
tion (5.277) for both the source and receiver. 

4. Consult Tables 5.1 and 5.2 and perform the indicated operation and obtain 
the limits of integration ^ 1 , 2 - 

5. Depending on type of partial haif-rim (whether the ends or the center of 
the half-rim is missing), evaluate the appropriate integrals and obtain the 
correction terms. 

5.6 The Diffraction-Point Method Solution with Corrections for Rim- 
Caustic Regions 

The corrections for the rim-caustic regions are provided by the transition 
functions in the diffraction coefficient. The transition function used depends on the 
rim-caustic region considered. For rim-caustic regions associated with the specular 
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Table 5.1: Table for evaluation of partial half-rim. 


Region 

Point 

Partial Rim 

r Operation 

|R* 

\Rr 




Lq4 


^cnd 

ID 




n 


n 

n 

A 

ViH 

0 

0 

lO 


D 

D 


B 

mm 


B 

1 

1 

0 

0 

0 

□ 

El 

n 

C 

n 


0 

KH 


1 

O 

Q 

El 

B 


D 

Q 


n 

0 

0 

0 

0 

0 

0 

H 

n 

E 

II 

n 

B 

0 

0 

1 

0 

0 

0 

B 

D 


O 

B 

mm 

0 

1 

0 

0 

0 

0 

B 

H 

A 


lEfl 

0 

mm 

mm 

■■ 

0 

Q 

n 


H 

B 

0 

0 

B 

B 

B 

■I 

1 

11 

H 

0 

B 

C 

n 

0 

B 

0 

0 

H 

0 

n 

0 

1 

B 

C 

D 

0 

0 

mm 

mm 

0 

0 

0 

n 

0 

B 

D 

n 

0 

0 

B 

B 

0 

0 


B 

D 

B 

E 

0 

■1 

1 

0 


1 

□ 


H 

9 

B 

E 

0 

n 

0 

mm 

n 

0 

□ 


n 

0 

B 

F 

VlH 

1 


B 

1 

0 

D 


0 

0 

C 

A 

mm 

0 

0 

0 

0 

mm 

O 


1 

0 

C 

B 

D 

0 

0 


0 

H 

0 


1 

D 

C 
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Table 5.2: Table for evaluation of partial half-rim. 
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or forward-scatter direction of the endcap, we use the function T^. For the other 
regions, the functions Tgi and Tg2 are used. The new diffraction coefficient with 
corrections for the rim-caustic region is 


D{Qe) = -$o$oD,{Qc) - 


where 


D,(Q>) = 


DM,) = 


— c 

2ny/2nk sin/9o 
-e~-^T 

2ny/2nkam0o 


{AT,i{U) 

(^Te2(t^) 


- BT,2{U)) 


+ btmu)) 


The transition functions are 

T,l(U) = ^(/cl(U) - 

T.2(t/) = \UciW) + 

U{Qe) = kaehe-{f + f') 


if (no -r') > 0 or (ho • f) > 0 then 


T.i(t/) = T,2(U) = Te(U) 


T.(t') = -/com 

7T 



f) 


where 


/co(U) 



i{u) + jJi{u)] -Ca(u) 


(5.291) 

(5.292) 

(5.293) 

(5.294) 

(5.295) 

(5.296) 

(5.297) 

(5.298) 

(5.299) 

(5.300) 
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(S.301) 


UW) 


-Cb(U) 


MV) 



+ J 




V 


) 


-Cc{U) 


(5.302) 


and Cji{U)j Cg{U), and Cc{U) are the correction factors for a partial half-rim. 
Outside the rim-caustic regions the transition functions are unity. 

The diffraction-point method solution with the rim-caustic correction, except 
for the new diffraction coefficient, has the original form which is given by 


« ^Lal = Ef + Ei + E{ + Ei (5.303) 

where 

Et = ^(«l) •5(g,)v//(<?, (5.304) 
^2 = ^{Ql) ■ S(«2)\//^(e2)«'‘"“ '^- (5.305) 

^3 = E^iQi) ■ B(«3)\/c''(93)<’‘''’ '^ (5.306) 

Ei = f*(Q4) .^(Q4)yy^eJ^4 . (5.307) 

Figure 40 shows how the calculated pattern, using the diffraction-point method 
with rim-caustic corrections, compares with the original diffraction-point method 
solution. 

The curved-side specular-region correction of the previous chapter can be 
introduced into this solution by multiplying the appropriate reflection-boundary 
term of the diffraction coefficient by the transition function (T). This is shown in 
the following diffraction coefficient with both the rim-caustic and curved-surface 
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specular-region corrections: 


where 

= cot ^ + cot ~ ' L~ ^ ^ ) (5.310) 

= cot (— T(±yXa,) +cot 

(5.311) 

The corrected diffraction-point method solution is 

^’»^<‘'olal=^l+®2+^3+^4 (5.312) 


where 

E{ = £•(<?,) (5.313) 

^2 = £'(<?2) ■B'(g2)v'/(«2)tl*'’*'^ (5.314) 

Bi = BHQ$) ■ D* (<?3) (5.315) 

= £’(04) • 5’ (Qt)\! ■ (5.316) 

The “corrected diffraction-point method solution'’ is the first-order UTD edge- 
diffraction solution with uniform corrections for the rim-caustic regions and for 
the curved-surface specular region. In Chapter 6, we will show how the corrected 
diffraction-point method solution compares with experimental backscatter mea- 
surements. And in Chapter 7, we will compare it against the moment method 
solution. 
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CHAPTER VI 


Comparisons with Experimental Measurements 


In order to verify our theoretical results and to provide confidence in our 
computer-generated patterns, a comparison with experimental measurements will 
be made. The measurements were performed in the Ohio State University Elec- 
troScience Laboratory compact radar range. 

Figure 41 shows the circular cone frustum that was used. It is constructed 
out of rolled brass sheet for the surface and brass plates for the endcaps which 
were soldered together. All seams were covered with metal tape. During the 
measurement, the cone frustum is held by a styrofoam column on top of a low RCS 
pedestal. The contributions from the mounting and background can be subtracted 
out of the measurements. 

Two types of measurements were performed on the cone frustum. One is 
a backscatter measurement at a fixed frequency as the frustum is rotated 360 
degreees — this generates a scattering pattern in the plane of rotation. The second 
type of measurement generates an impulse response. Backscatter mesisurements 
for a fixed aspect angle are taken, as the frequency is swept over a broad range. 
Performing an inverse Fourier transform on the data produces a band-limited ver- 
sion of the impulse response. 

The principal plane backscatter measurements are shown in Figures 42 and 43 
which are compared with the calculated UTD results. The measurements were 
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Figure 41: Cone frustum used for measurements, / = 4.00 GHz. 
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Figure 42: Cone frustum backscatter pattern for horizontal polarization at 4.00 
GHz obtained using UTD (solid line) and measurement (dashed line). 
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Figure 43: Cone frustum backscatter pattern for vertical polwization at 4.00 
GHz obtained using UTD (solid line) and measurement (dashed line). 
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performed at a frequency of 4 GHz for both horizontal and vertical polariztions. 
Symmetry in the measured pattern gives a good indication of its reliability. There 
is excellent correlation at most of the significant peaks (within 1 dB), and the 
differences can be explained. 

Higher-order diffractions are probably the best explaination for the differences 
between the measured and calculated patterns. The effects of double edge diffrac- 
tion are usually strongest at the surface boundaries and for the horizontal (hard) 
polarization case. (The vertical polarized wave tends to get shorted out along the 
conducting surface.) Looking at the pattern for horizontal polarization, we note 
the discontinuities and the differences in the calculated result at 90 and 270 de- 
grees. This indicates that multiple diffractions across the endcaps are significant 
in this region. The discrepancies at 15 and 345 degrees are also probably due to a 
higher order diffraction along the curved surface. 

Higher-order terms in the UTD solution may also resolve the difference at 0 
degrees a::ial incidence. We notice that there is about a 5 dB difference. At 0 
degrees, we are facing the frustum nose-on, and it is a very frequency sensitive 
region due to the interference between the diffracted fields of bax:k and front rims. 
Let us see what happens when the frequency of the calculated pattern is shifted 
slightly from 4.00 GHz to 4.05 GHz (see Figures 44 and 45). An increase of 5 dB 
occurs at 0 degrees with very little change in the pattern elsewhere. The reason for 
such a dramatic change will be clearer when we look at the backscatter frequency 
scan plot. 

Figure 46 shows the calculated UTD frequency scan plot for 0 degrees axial 
incidence. (Polarization does not matter due to the symmetry.) The oscillation 
observed in the plot is due to the interaction between the contributions from the 
front endcap and the back rim. As the frequency decreases, the variation in the 
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Fieure 44: Cone frustum backscatter pattern for horizontal polarization at 4.05 

GHz obtained using UTD. 
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Figure 45: Cone frustum backscatter pattern for vertical polarization at 4.05 

GHz obtained using UTD. 
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Figure 46: Calculated backscatter plot using UTD for a scan in frequency from 2 

to 18 GHz at 0** axial incidence. 
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magnitude becomes greater due to the increasing diffracted field from the back rim 
over the specular reflection from the front endcap. Note that neau* 4 GHz we have 
a sharp null. This explains why we have a large change in magnitude with only a 
slight change in frequency. Higher-order terms, which become more significant at 
lower frequncies, might be responsible for slight phase differences that would shift 
the actual backscatter frequency scan plot. 

Next, let us look at the backscatter impulse responses at 0 degrees, or axial 
incidence. The frequency was varied from 6 to 18 GHz. The calculated and 
measured plots are shown in Figures 47 and 48. An inverse Fourier transform is 
performed and the resulting band-limited impulse reponses are shown in Figure 49. 
The phase references of measured and calculated patterns are not the same — which 
results in a shift in time. The impulse response from the measured data shows a 
dominant contribution from the from endcap followed by a contribution from the 
back endcap rim. This is followed by a lesser contribution which is a double 
diffraction across the back endcap. The calculated plot compares well, except of 
course, there are no double diffractions. 
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Figure 47: Measured backscatter plot (UTD dashed line) for a scan in frequency 

from 6 to 18 GHz at 0** axial incidence. 
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Figure 48: Calculated backscatter plot using UTD for a scan in frequency from 6 

to 18 GHz at 0** axial incidence. 
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Figure 49: Impulse response using measured (top) and calculated UTD (bottom) 
backscatter data from 6 to 18 GHz for 0® axial incidence. 
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CHAPTER VII 


Comparisons with the Moment Method Solution 


In the last chapter comparisons were made with backscatter measurements. 
The good agreement gives us greater confidence that the theory and computer 
code are working as it should. However, the theory has not yet been tested for 
bistatic angles. 

In this chapter, the calculated UTD results will be compared with results 
obtained by using the moment method. (We have a great deal of confidence in 
the moment method solution which is regarded as an “exact” solution for these 
calculations.) We make our comparisons with the moment method solution because 
accurate bistatic measurements are often difficult to obtain. This is especially true 
for large bistatic angles where there might be significant coupling of the incident 
field into the receiver. 

For the cone frustum, the moment method solution presented in reference 
[16] will be used. It can determine the plane-wave scattering for any rotationally 
symmetric conducting body. Using this solution, C.W. Chuang [17] has developed 
a computer code that can generate bistatic scattering patterns for fixed angles of 
incidence. 

The fixed incidence angles of = 0, 15, and 30 degrees are chosen for the 
cone frustum. Larger angles are possible, but due to the longer computer run 
times required, they are avoided. The frustum dimensions will be the same as that 
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used for the measurements (see Figure 41). Again, due to the amount of computer 
time, larger dimensions were not chosen. 

Now let’s look at the calculated bistatic principal plane patterns. Figures 50 
through 55 show the comparison between the results obtained using UTD and 
the moment method solutions for both horizontal and vertical polarizations. 

For most of the patterns, agreement is within 3 dB or better. In the specular 
and forward'Scatter regions, the match is within 1 dB or better. However, near 
backscatter for the 0 degree axial incidence pattern, we have a significant difference. 
This is a rim-caustic region which involves both front and back rims. We note that 
this the same region where we had previous differences with the measurements. As 
discussed in the last chapter, higher-order diffractions may need to be considered 
to resolve this difference. 

Higher-order diffractions probably account for the other differences between 
the UTD and MM patterns. Double diffractions are significant along the surfaces 
boundaries at 30, 150, 90, 210, 270, and 330 degrees due to interactions across the 
endcaps and between front and back rims. The creeping- wave and creeping- wave 
interactions with the edges are probably significant contributors in the forward- 
scatter direction. But it is dangerous to speculate without further investigation. 
Overall, the agreement is reasonable for engineering purposes. At higher frequen- 
cies, the accuracy of the first-order UTD solution will get better, due to the greater 
attenuation of surface waves that cause higher-order effects. 

Finally, one might wonder why we need a UTD solution when a more accurate 
moment method solution is available. The relative calculation times might be 
worth noting: the UTD result took 0.3 minutes of CPU time to generate the 
pattern; the moment method result took 3 hours (on a VAX 11/780). The faster 
calculation time seems to be a reasonable trade-off for accuracy. 
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Figure 50: Bistatic principal plane pattern for fixed incidence at d — 0 degrees 
and horizontal polarization, obtained using the moment method (dashed line) 

and UTD (solid line). 
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Figure 51: Bistatic principal plane pattern for fixed incidence at 9 = 0 degrees 
and vertical polarization, obtained using the moment method (dashed line] and 

UTD (solid line). 
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Figure 52: Bistatic principal plane pattern for fixed incidence at • = 15 degrees 
and horisontal polarization, obtained using the moment method (dashed linej 

and UTD (solid line). 
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Figure 53: Bistatic principal plane pattern for fixed incidence at = 15 degrees 
and vertical polarization, obtained using the moment method (dashed line) and 

UTD (solid line). 
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Figure 54: Bistatic principal plane pattern for fixed incidence at 0 = 30 degrees 
and horizontal polarization, obtained using the moment method (dashed line) 

and UTD (solid line). 
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Figure 55: Bistatic principal plane pattern for fixed incidence at = 30 degrees 
and vertical polarization, obtained using the moment method (dashed line) and 

UTD (solid line). 
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CHAPTER Vra 


Summary and Conclusion 

The objective of this work was to solve for the high-frequency scattered fields 
from a cone frustum. The desired solution needed to be practical — suitable for 
engineering applications. Thus, we chose a solution based on the methods of the 
geometric theory of diffraction — for its accuracy, for the physical insight that ray 
methods provide, and for its ease of application. The original GTD solution for the 
frustum, which we called the “diffraction-point method solution”, provided fairly 
accurate results, but it failed in the specular region for the curved surface and in 
the rim-caustic regions. 

We wanted to correct the original GTD solution in these regions, but we 
also wanted to avoid direct switching to alternate solution forms, such as the 
physical optics solution. The diffraction-point method solution has a simple and 
elegant form, and the virtues of a simple solution can be appreciated in applications 
dealing with complicated systems and computer codes. Therefore, to retain the 
original form, we developed correcting functions, called “transition functions” , that 
modified the terms in the diffraction coefficients. The transition functions provide 
a uniform solution — in other words, they are modified so that they can be used 
in all regions. And as demonstrated in the last two chapters, they provide results 
that compare well with measurements and moment method results. 

We have accomplished our objective, but many items still need further inves- 
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tigation. The following are a few questions and items that need more study: 


• Rim-caustic transition functions. A single more general function is needed 
that works well for all regions, for elliptical as well as circular rims. 

• Quasi-near-zone scattering. Can we model the behavior in the reflection 
caustic regions from far-zone to near-zone? Would it be possible to develop 
another transition function to do this? 

• Scattering outside the principal plane. Do other mechanisms, such as creeping- 
wave interactions with the edge, become important? 

• Higher-order diffraction effects. The development of solutions that consider 
muliply diffracted rays may extend the UTD solution to lower frequncies. 

• Scattering from structures using cone frustums. Some of the applications, 
such as the scattering from multiple cone-frustum structures in combina- 
tion with other geometries can be investigated. What shadowing and ray 
interactions between structures need to be considered? 

The author hopes that some of these questions will be answered in future work, 
and hopefully this report will serve as a basis for that research. 
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APPENDIX A 


Confluent Hypergeometric Functions 

i 

The confluent hypergeometric functions (also called Hummer functions) are 
solutions to the confluent hypergeometric equation, 

zy"(z) + (6 - z)y'{z) - ay = 0. (i4.317) 

I 

The complete solution is given by, 

y{z) = AM{a,b,z) + BV{ayb,z) (A.318) 

where A and B are are arbitrary constants and a ^ — n, 6 ^ —m (m, n are positive 
integers). 

j The confluent hypergeometric functions M and U have integral, series, and 

I 

large argument forms. References for this section are [13] and |18|. The integral 
representation is given by 

i 

Re{b) > Re{a) > 0 

i 

j 

' U(a,4,s) = ^ + (A.320) 

r(a) Jo 

Re(z)>0,Re{a)>0 
The series form of the functions are 
M(a,6,r) = 

' 148 


I 


(o)n = o(o + l)(o + 2)...(a + »-l) 
(o)o = 1 


U(a,6,z) 


K ( M(a, 6, z) 

sin7r6 tr(l + a-6)r(6) 

^ r(a)r (2 - 6) / 


(A.322) 


The large argument form, for \z\ large, (a and b fixed) is 


M(a,6, 2 ) = 


r(t) 

r(i - <.) 




A-1 

E 

n=0 


(o)n(l o ~ b)n , 


(-»)-" + o(i»r*) 


n! 


] 




r(o) 



(6-.)„(l-»)„ ^_„ ^ p i 

ft! 


1 


(A.323) 


where the sign is (+) if — jTT < argz < jx, or (— ) if — jtt < argz < 


U(a,6,z) = z 


-a 


n=0 


(A.324) 


where — §?r < arga < j?r. 

The two confluent hypergeometric functions of interest in this work are M(l, 2.5, -z) 
and U(l, 2.5, -z). For small values of the argument, the series form should be used. 

The series form is given by 


M(l,2.5,-z) 



(A.325) 
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(A.326) 


U{1,2.5,-*) 


-?M(l,2.5,-z) + -^ 

® 2(-*)t 


For large arguments, \z\ large, 


* .5)n 


Z2 ^ n=0 

1 

U(1,2.5,-2) = — 5^{-.5)„z- 


n=0 


Usually three terms, 5 = 3, are sufficient. 


(A.327) 

(A.328) 
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APPENDIX B 


Bessel and Struve Functions 

In Chapter 5, we use the zero and first order Bessel and Struve functions. 
These functions can be calculated by using a power series for small values of the 
argument, or by using the asymptotic form for large values of the argument. The 
references used are [18,19]. 

The power series expansions are given by the following: 


Jn(x) 

(-1)*(x/2)2*+" 
^0 *!r(* + l + n) 

(B.329) 

M^) 

ki, (*')* 

(B.330) 

Jl(x) 


(B.331) 

Hn(x) 

(-l)‘(*/2)“ 

^ jPo + 1.5)r(fc + n + 1.5) 

(B.332) 

Ho(x) 

- ^ ,\ib+l_2ifc-l f 2*^’] 

- ^ * [wi] 

(B.333) 

Hi(x) 

- -y-l) X (2t + I)[,j(^^,„,J . 

(B.334) 
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In terms of the confluent hypergeometric function, the Bessel function is given by; 


Jn(x) = ■ 

The following are the asymptotic expansions for large |a:| : 


, , X r , “1“2 . 01030507 \ 

~ ) 

. , fOLl O1Q3O5 , >^1 


where 




O1O3O5O7 

4! 


_...) 


, /ai 010305 , 'll 
+ cosV.(^-— 55 — + •••]] 


, nir ir 

xb = X 

^ 2 4 


Ofc = 


4n^ — 
8x 


and 


Ho(x) ~ 


1 1 l2 . 32 l2 . 32 • 5 ^ 


3 + 5 

X X^ X* 


7T 

+l'o(*) 


Hi(i) ~ - 


1 + “7? “ 


12.3 l2 . 32 . 5 
~TT~ ^6 


+ri(x) . 


(B.335) 


(B.336) 

(B.337) 

(B.338) 

(B.339) 

(B.340) 

(B.341) 
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APPENDIX C 


Computer Programs 

This appendix contains some of computer subroutines and functions that were 
used to calculate the confluent hypergeometric functions and the integrals used for 
the partial rim corrections. The programs are written in FORTRAN 77. They are 
based on the material found in Chapter 5 and Appendix A. 
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c* ************************************************ 
C This function is th« confluent 

C hyper geometric function of 

C U(1.,2.5,Z) in series form. 

C Used with small values of the 

C argument (magnitude of Z < 5 ) . WE 8/9/85 

C 

COMPLEX FUNCTION UPS(Z) 

C 

C Z > function argument, complex 

C 

COMPLEX CJ,MFS,Z,Z1,Z2 

C 

CJ-(0.,1.) 

PI-3.14159265 

C 

Zl->(2./3. )*MFS(1. ,2.5,Z) 
Z2-.5*SQRT(PI)*CEXP(Z)/(Z**1.5) 

C Z2 can become infinite if Z-0. 

UFS-Z1-«‘Z2 

RETURN 

END 

C* ************************************************ 

C This function is the confluent 

C hypergeometric function of 

C U(1.,2.5,Z) in asymptotic form. 

C Used with large values of 

C the argument (magnitude of Z >> 1 ). 

C 

COMPLEX FUNCTION UFL(Z) 

C 

C Z - function argument, complex 

C 

COMPLEX S,Z 

C 

S-(0. ,0. ) 

c 

C J is the number of terms. 

J-3 

C 

DO 10 N-0,J-1 
X-N 

S-(PH(-.5,N)/(-Z)**X)+S 
10 CONTINUE 

UFL-S/Z 



RETURN 

END 
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C************************************************ 

C This function is used with functions 4JFL. 

C (It is often represented by the 

C Pochhamner syabol.) 

c 

FUNCTION PH(A,N) 

P-1. 

IF(N.EQ.O)TBEN 

PH-1. 

ELSE IF(N.EQ.1)THEN 
PH-A 

ELSE IF(N.GE.2)THEN 
DO 10 K«0,N-1 
P-(A+K)*P 

10 CONTINUE 

PH-P 
ENDIF 



RETURN 

END 

C* ***************************************** ***♦*< 

C This function is the confluent 

C hyper geometric function of H(A,B,Z) in 

C series form. Used with small values of 

C the argument (magnitude of Z < 5 ) . 

C WE 8/9/85 

C 

COMPLEX FUNCTION MFS(A,B,Z) 



C A,B - parameters of the hypergeometric 

C function (contents) 

C Z - function argument, complex 



DIMENSION CN(0:40) 

COMPLEX SS,Z 



SS-(0.,0.) 



X-CABS(Z) 



C K - number of terms used in the series 

C 

K-3.*X+11. 

KMAX-40. 

1F(K.GT.KMAX)K-KHAX 
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c 

C Small argument. 

C 

ir(X.LT..01)NFS-(l.,0.) 

c 

C Series form. 

C 

IF(X.GE. .ODTHEN 

CALL CS(A,B.K.CN) 
J-K-1 

DO 10 N-l.J 

M-J-(N-l) 

SS-Z*(CN(H)-fSS) 

10 CONTINUE 

SS-1.4-SS 
HFS-SS 
ENDIF 

C 

RETURN 

END 


C This subroutine calculates the constants 

C used in the series form of the 

C hypergeometric function M(a,b,jX). 

C 

SUBROUTINE CS(A,B,K,CN) 

C 

C A.B > parameters of the hypergeometric 

C function 

C K ■ number of terms in the series 

C CN array of constants to be returned 

C 

DIMENSION CN(0:40) 

C 


CN(0)-1. 


DO 10 N-1,40 
X-N 


IF(N.LE.K)THEN 

CN { N ) -CN ( N-1 ) * ( A+X-1 ) /( ( B+X-1 ) *X ) 
ELSE 

CN(N)-0. 

ENDIF 

10 CONTINUE 

RETURN 
END 
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C*************************************«******«* 

This function is the confluent 
hypergeometric function of M(1.5r2.5«Z) 
in asymptotic form. Used with large 
values of the argument 
(magnitude of Z >> 1 ). WE 85 

COMPLEX FUNCTION MFL(Z) 


Z • function argument, complex 


COMPLEX C1,C2,CJ,S,Z 


Pl-3. 14159265 
CJ-(0. ,1. ) 
S-(0. ,0. ) 


ZI-AZMAG(Z) 

ZM-CABS(Z) 

IF(ZI.EQ.-ZM)C1-.75*SQRT(FI) 
1 *CEXP(-CJ*1.5*PI)/(Z**1.5) 

IF(ZI.EQ.ZM)C1-.75*SQRT(PI) 

1 *CEXP(CJ*1.5*PI)/(Z**1.5) 


J is the number of terms. 
J«3 

C 

DO 10 N-0,J-1 
X-N 

S-(PH(-.5,N)/Z**X)+S 
10 CONTINUE 

C2-(1.5*CEXP(Z)/Z)*S 

MFL-C1+C2 

C 

RETURN 

END 
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c* ************************************************ 
This subroutine approxinately computes 
the integrals used in the transition 
functions from XI to X2, for small 
values of U. 

SUBROUTINE INTGRL ( XI , X2 , U , NINT , PX , ECODE ) 


Xl,X2 ■ limits of integration 
U - variable in the exponential 

(a real number with mag. < 8.) 

NZNT ■> 1 cos(x) * exp(j u cos(x)) 

2 ( cos(x)**2 ) * exp(j U cos(x)) 

3 ( sin(x)**2 ) * exp(j U cos(x)) 
FX - computed value of the integral 
ECODE ■ error code 

The subroutines INCOS and INSINCOS are 
used with this routine. 


COMPLEX CJ,FX,SUM 

DIMENSION FI(0:25),C(0:23) ,FACT(0:23) 

C 

DATA PI, TPI, DR/3. 14159265,6. 28318537, 

1 0.017453292/ 

DATA FACT/1 . , 1 . , 2 . , 6 . , 24 . , 120 . , 720 . , 5040 . , 

1 40320., 362880. ,3628800., 39916800., 

2 479001600. ,6227020800. ,8. 7178291E10, 

3 1.3076743E12,2.0922789E13,3.5568742E14, 

4 6.4023737E15,1.2164510E17,2.4329020E18, 

5 5.1090942E19,1.1240007E21,2.5852016E22/ 
CJ-(0. ,1. ) 

C 

XMIN-.OOl 

NNAX-23 

C 

IF(ABS(U) .GT.8. )THEN 
ECODE-4. 

GO TO 999 
ENDIF 

C 

DO 10 K-0,23 

C(K)-U**K/FACT(K) 

IF(ABS(C(K) ) .LT.XMIN)THEN 
NMAX-K 
GO TO 20 
ENDIF 

10 CONTINUE 
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20 CONTINUE 

sun-(o. ,0. ) 

XF(NXNT.EQ.1)THEN 

CALL XNCOS(Xl,X2,NHAX-fl,FX) 

DO 30 K-0,NHAX 

sun-sun-t- ( cj* *k ) *c ( k ) *fx ( K-t>i ) 

30 CONTINUE 

ELSE XF(NXNT.EQ.2)TBEN 

CALL XNCOS(Xl,X2,NnAX-t>2,FX) 

DO 40 K-0,NKAX 

SUM-SUM+ ( C J**K ) *C ( K ) *FX ( K+2 ) 

40 CONTINUE 

ELSE IF(NINT.EQ.3)THEN 

CALL INSINCOS(Xl,X2,NHAX,FI) 

DO 50 K-0,NMAX 
SUM-SUM+(CJ**K)*C(K)*FI(K) 

50 CONTINUE 

END IF 

C 

FX-SUM 

999 CONTINUE 

RETURN 
END 

C**** *************** ****************************** 

C This subroutine computes the integral of 

C the cosine to the power N. It computes up 

C to a specified value of N and returns 

C the result in an array where the array 

C subscript corresponds to the value of N. 

C 

SUBROUTINE INC0S(Xl ,X2 ,NHAX,F1 ) 



C XI, X2 - limits of integration (radians) 

C NMAX - maximum value of power N 

C (an integer) 

C FI ■ the array of computed integrals 



DIMENSION FI ( 0:25 ) ,F1( 0:25) ,F2( 0:25) 



F1(0)-X1 

F2(0)-X2 

FI(0)-X2-X1 

FKl)-SIN(Xl) 

F2(1)-SIN(X2) 

FI(1)-F2(1)-F1(1) 

F1(2)-SIN(2.*X1)/4.+X1/2. 

F2(2)-SIN(2.*X2)/4.+X2/2. 

FI(2)-F2(2)-F1(2) 
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XF(NHAX.GT.2)THEN 
DO 10 K-3,NMAX 

A-K 

ri(K)-(SIN(Xl)*(COS(Xl)**(K-l))/A)+ 

1 ( (A-1. )*Fl(K-2)/A) 

F2(K)-(SIN{X2)*(COS(X2)**(K-l) )/A)+ 

1 ( (A-1. )*F2(K-2)/A) 

FI(K)-r2(K)-rl(K) 

10 CONTINUE 

ENDIF 
RETURN 
END 

C* ******** **************************************** 

This subroutine computes the integral of 
the cosine to the power N multiplied by 
the sine squared. It computes up to a 
specified value of N and returns the 
result in an array where the array 
subscript corresponds to the value of N. 

SUBROUTINE ZNSINCOS (Xl ,X2 ,NHAX, FI ) 


XI, X2 - limits of integration (radians) 
NMAX - maximum value of power N 
(an integer) 

FI “ the array of computed integrals 


DIMENSION FI(0:25) ,F1( 0:25) ,F2 (0:25) , 
1 FCOS(0:25) 

C 

CALL 1 NCOS ( XI, X2, NMAX, FCOS) 

C 

DO 10 R-0,NHAX 
A-K 

F1(K)— SIN(X1)*C0S(X1)**(K+1) 

F2(K)— SIN(X2)*COS(X2)**(K+l) 

FI (K)-( (F2(K)-F1(K) )+FCOS(K) )/(A+2. ) 
10 CONTINUE 

C 

RETURN 

END 
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